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Our topic concerns a long standing puzzle: the energy of gravitating systems. 
More precisely we want to consider, for gravitating systems, how to best describe 
energy-momentum and angular momentum/center-of-mass momentum (CoMM). 

It is known that these quantities cannot be given by a local density. The modern 
understanding is that (i) they are quasi-local (associated with a closed 2-surface), 

(ii) they have no unique formula, (iii) they have no reference frame independent 
description. In the first part of this work we review some early history, much 
of it not so well known, on the subject of gravitational energy in Einstein’s gen¬ 
eral relativity (GR), noting especially Noether’s contribution. In the second part 
we review (including some new results) much of our covariant Hamiltonian for¬ 
malism and apply it to Poincare gauge theories of gravity (PG), with GR as 
a special case. The key point is that the Hamiltonian boundary term has two 
roles, it determines the quasi-local quantities, and, furthermore it determines the 
boundary conditions for the dynamical variables. Energy-momentum and angular 
momentum/CoMM are associated with the geometric symmetries under Poincare 
transformations. They are best described in a local Poincare gauge theory. The 
type of spacetime that naturally has this symmetry is Riemann-Cartan space- 
time, with a metric compatible connection having, in general, both curvature and 
torsion. Thus our expression for the energy-momentum of physical systems is 
obtained via our covariant Hamiltonian formulation applied to the PG. 


1. Introduction 

How to give a meaningful description of energy-momentum for gravitating systems 
(hence for all physical systems) has been an outstanding fundamental issue since 
Einstein began his search for his gravity theory, general relativity (GR). It is deeply 
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connected to the essential nature of not only geometric gravity but of all the fun¬ 
damental interactions—their inherent gauge nature. Noether’s paper that includes 
her two famous theorems relating global symmetries to conserved quantities and 
local gauge symmetries to differential identities was originally motivated by this 
very issue. She showed that gravitational energy has no proper local description. 
So investigators only found various expressions which were inherently non-tensorial 
(reference frame dependent pseudotensors). They have two inherent ambiguities: 
(i) there are many possible expressions, (ii) they are non-covariant— reference frame 
dependent. The modern view is that energy-momentum is quasi-local (associated 
with a closed 2-surface). Quasi-local proposals have analogous ambiguities. These 
ambiguities can be clarified by the Hamiltonian approach. From a first order La- 
grangian for quite general differential form fields, we have constructed a spacetime 
covariant Hamiltonian formalism, which incorporates the Noether conserved cur¬ 
rents and differential identities. The Hamiltonian that dynamically evolves a spatial 
region includes a boundary term. The explicit form of the boundary term depends 
on the boundary conditions and also on an appropriate reference choice. With 
a suitable vector field it gives expressions for the quasi-local quantities (energy- 
momentum, angular momentum/center-of-mass momentum) and also quasi-local 
energy flux. A geometric gauge theory perspective provides the most appropriate 
dynamical variables. The geometry is Riemann-Cartan, with, in general, both cur¬ 
vature and torsion. For the PG (GR is a special case) with general source and 
gauge fields, we identified a preferred Hamiltonian boundary expression along with 
a procedure for finding a “best matched” reference. With this one can obtain val¬ 
ues for the quasi-local energy-momentum and angular momentum/center-of-mass 
momentum |fj 

Our topic here concerns the localization of energy-momentum. The main aim 
of our research program has actually been to better understand the Hamiltonian 
for dynamic spacetime geometry, especially the role of the Hamiltonian boundary 
term. It turns out that this sheds much light on the issue of the localization of 
energyPS A number of different ideas will be fit together to give a good picture of 
this long standing puzzle. In addition to being mindful of Noether’s results, we will 
use a Hamiltonian approach combined with a local gauge theory view of dynamic 
spacetime geometry. 

This present work is largely just an application of Noether’s result. We will 
begin with some early history (much of it not so well known) regarding energy in 
the context of GR, especially Noether’s contribution. Next, we will show how in 
GR pseudotensors are connected with the Hamiltonian and introduce the quasi¬ 
local idea. Following this we make some brief remarks about gravity, geometry, 
connection and gauge. We then introduce and give a short review of our main tool: 
differential forms. We develop in some detail variational principles with differential 


a For an alternative to our Hamiltonian approach to energy-momentum and angular momentum 
for the PG see Kawai (1986).^ 
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form fields. With this we can give simple examples of applications of the two 
Noether theorems. We then introduce the first order formalism followed by our 
Hamiltonian formulation and the 3+1 split. The Hamiltonian boundary term and its 
important roles are discussed next. Asymptotic fall offs for the fields are noted. We 
explain why Riemann-Cartan geometry is appropriate for our purposes. Variational 
principles for dynamic spacetime geometry with quite general sources are developed, 
including the Noether conserved currents and differential identities. We present a 
first order and Hamiltonian formulation for the PG along with the Hamiltonian 
boundary term and identify our preferred Hamiltonian boundary term for these 
dynamic spacetime geometry theories. We also include a prescription for choosing 
the necessary reference values that are needed for the quasi-local energy-momentum 
and angular momentum expressions. 

2. Background 

As this present work approached its final form we received some very good news: 
all of the volumes of the Einstein papers published so fai0 both the originals and 
the English translations 3 ®—are now freely available online.^ An examination of 
Einstein’s dozens of papers on gravity during the period 1913—1918, as well as his 
extensive correspondence with his contemporaries on the topic of gravity, shows 
that most of them include a significant consideration of the topic of gravitational 
energy. 

2.1. Some brief early history 

We have only begun to look into the historical development of the modern ideas 
regarding gravitational energy; the topic merits much deeper study. Here we can 
only give a brief report, relying on the Einstein papers as well as some of the many 
good historical investigations available, in particular regarding energy-momentum 
conservation! 13 1 18 !^ 1 

We will rely on the Hamiltonian formalism applied to dynamical variables that 
are related to a local gauge theory of spacetime symmetry approach. (EarmaiP® 
has given a very interesting discussion on how the Hamiltonian approach connects 
with the gauge theory perspective.) 

It seems not so well known that gravitational energy, or more precisely the proper 
description of the energy of gravitating systems (i.e., all real physical systems), has 
played a large role in the development of 20th century physics. 

In the years 1912-1915 Einstein, when he was searching for satisfactory field 
equations, used a form of the equations that explicitly included an energy- 
momentum density for the gravitational field and were designed to satisfy the 
principle of energy-momentum conservation @ Thus an expression for the Einstein 

b At present up to 1923. 

c See Refs. I61H63I94I98I126I for discussions of how Einstein found his field equationsPU 
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energy-momentum pseudotensor already existed even before he found the correct 
field equations. It should be appreciated that general covariance brought with it 
features that had never before been encountered in any theory. (Indeed there is still 
controversy up to the present day!— I— I 96 ! 125 ! ) For a couple of years Einstein very 
much doubted that a generally covariant theory could be foundjj he proposed that 
energy conservation would select the preferred physical coordinate frame. Initially 
Hilbert followed Einstein in this belief (see the proofs of his first note in Ref. [108] '). 

Although Einstein began using variational principles in 1914231 this was not his 
path to the field equations. Hilbert was the first to identify a generally covariant 
Lagrangiar@ (proportional to the Riemannian scalar curvature). He also constructed 
(in a complicated way that was not easy to understand) his “conserved energy 
vector”, a vector with vanishing divergence associated with the general coordinate 
invariance (i.e., diffeomorphism invariance) of his Lagrangian. 

Einstein’s energy-momentum pseudotensor was criticizecP 7 for giving “unphys¬ 
ical” values (Schrodinger 11 ^ noted that one could choose the coordinates to give a 
vanishing value outside a fluid sphere, and Bauei® noted that one could choose the 
coordinates to give a vanishing energy value for empty Minkowski space). 

Lorentz, Levi-Civita and Klein argued that the Einstein curvature tensor G 
was the only proper gravitational energy-momentum density; hence one should 
regard the Einstein equation in the form 


- Gfj, v + T^ u — 0 ( 1 ) 

K 

as describing the vanishing sum of gravitational and material energy-momentum. 
(This idea has been advanced more recently by Cooperstockf^l) In our modern 
perspective for GR their idea is quite correct—but a density is not the whole story. 
There is more to energy-momentum than just a density. 


2.2. From Einstein’s correspondence 

Here are some excerpts from Einstein’s correspondence concerning the Einstein 
pseudotensor, Hilbert’s energy vector and Noether’s contribution! 13 ! 14 ! 16 ! 75 ! 92 ! 1 -^ 
They reveal the difficulties and the extent of understanding these people had at 
that time. All these are quoted from the Einstein paper JMEU Vol. 8. 

“Highly esteemed Colleague, ... I am sitting over your relativ¬ 
ity paper, ..., and am honestly toiling over it. I do admire your 
method, as far as I have understood it. But at certain points I 
cannot progress and therefore ask that you assist me with brief 
instructions. ... I still do not grasp the energy principle at all, not 
even as a statement.” (Doc. 221 to Hilbert 25 May 1916) 

d One reason was his famous “hole” argument ! 81 ! 98 ! 

e Einstein and Hilbert had quite different agendas! 112 ! 115 ! 136 ! 141 ! Hilbert in his Foundation of 
Physics papers, based on the work of Einstein and Mie, was using his axio matic method with 
the objective of finding a unified field theory of gravity and electromagnetism! 15 l 28 l 58 l lu8 l lu9 l 114 l 
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“Your explanation of equation (6) of your paper was charming. 

Why do you make it so hard for poor mortals by withholding the 
technique behind your ideas? ... In your paper everything is un¬ 
derstandable to me now except for the energy theorem. Please do 
not be angry with me that I ask you about this again. ... How is 
this cleared up? It would suffice, of course, if you would charge 
Miss Noether with explaining this to me.” 

(Doc. 223 to Hilbert 30 May 1916) 

“My tg’s are being rejected by everyone as unkosher.” 

(Doc. 503 to Hilbert 12 April 1918) 

.. Only (24) is an identity ... The relations here are exactly anal¬ 
ogous to those for nonrelativistic theories.” 

(Doc. 480 to Klein 13 March 1918) 

“I have succeeded in discovering the organic formation law for 
Hilbert’s energy vector” (Doc. 588 from Klein 15 July 1918) 

“The only thing I was unable to grasp in your paper is the conclu¬ 
sion at the top of page 8 that e a was a vector.” 

(Doc. 638 to Klein 22 Oct 1918) 

“Thank you very much for the transparent proof, which I under¬ 
stood completely.” (Doc. 646 to Klein 8 Nov 1918). 

“... Meanwhile, with Miss Noether’s help, I understand that the 
proof for the vector character of e a from “higher principles” as I 
had sought was already given by Hilbert on pp. 6, 7 of his first note, 
although in a version that does not draw attention to the essential 
point.” (Doc. 650 from Klein 10 Nov 1918) 

Briefly, after a couple of years Klein clarified Hilbert’s energy-momentum “vec¬ 
tor” ; he related it to Einstein’s pseudotensor, but (as we will discuss in more detail 
shortly) disagreed with Einstein’s physical interpretation of divergenceless expres- 
sions|j Enlisted by Hilbert and Klein, it was Emmy Noether who solved the primary 
puzzle regarding gravitational energy. 

2.3. Noether’s contribution 

If one had to describe 20th century physics in one word a good choice would be 
symmetry. Most of the new theoretical physics ideas involved symmetry. Essen- 

f For these investigators these things were not as easy as they are for us t oday; in particular the 
Bianchi identity and its contracted version were not known to these people p 8 ! 113 ! so they had, in 
effect, to rediscover that identity from effectively the diffeomorphism invariance of the Lagrangian. 
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tially they can be seen as applications of Noether’s theorems. Briefly, Noether’s 
first theorem associates conserved quantities with global symmetries, and Noether’s 
second theorem concerns local symmetries: it is the mathematical foundation of 
the modern gauge theories. Unfortunately her worlf] was largely overlooked for 
about 50 years Why did Noether make her investigation? To clarify the issue of 
gravitational energy. 

Klein was looking into Einstein’s theory and the relationship between Einstein’s 
pseudotensor and Hilbert’s energy vector. Some of the correspondence between 
Hilbert and Klein was published in a paper.^l We quote some excerpts^] 


Klein wrote 


You know that Miss Noether advises me continually regarding 
my work, and that in fact it is only thanks to her that I have 
understood these questions. When I was speaking recently to Miss 
Noether about my result concerning your energy vector, she was 
able to inform me that she had derived the same result on the basis 
of developments of your note (and thus not from the simplified 
calculations of my section 4) more than a year ago, and that she 
had then put all of that in a manuscript (which I was subsequently 
able to read). She simply did not set it out as forcefully as I recently 
did at the Mathematical Society (22 January [1918]). 

Hilbert responded 

I fully agree in fact with your statements on the energy theo¬ 
rems: Emmy Noether, on whom I have called for assistance more 
than a year ago to clarify this type of analytical questions concern¬ 
ing my energy theorem, found at that time that the energy compo¬ 
nents that I had proposed—as well as those of Einstein—could be 
formally transformed, using the Lagrange differential equations (4) 
and (5) of my first note, into expressions whose divergence vanishes 
identically, that is to say, without using the Lagrange equations (4) 
and (5). 

Also 


Indeed I believe that in the case of general relativity, i.e., in 
the case of the general invariance of the Hamiltonian function, the 
energy equations which in your opinion correspond to the energy 
equations of the theory of orthogonal invariance do not exist at all; 

I can even call this fact a characteristic of the general theory of 
relativity. 

s For discussions see Refs. 113114I16I75I921 1111 . 

h We do not know of any English translation of Klein’s papers; the translations of the following 
excerpts are quoted from Ref. I75| p 66. 
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What Hilbert here calls the Hamiltonian function we now refer to as the La- 
grangian. Noether wrote her 1918 paper to clarify the situation. 

2.4. Noether’s result 

Many have heard of Noether’s theorems, but the full scope of what she actually 
did is not so generally well known. So we take this opportunity to quote (from the 
highly recommended book of Kosmann-SchwartzbaclP^) in full her key results. (It 
should be mentioned that her Lagrangians were quite general, they could depend 
on any finite number of derivatives.) 

Theorem I. If the integral I is invariant under a [finite contin¬ 
uous group with p parameters] G p , then there are p linearly in¬ 
dependent combinations among the Lagrangian expressions which 
become divergences—and conversely, that implies the invariance of 
I under a [group] G p . The theorem remains valid in the limiting 
case of an infinite number of parameters. 

Theorem II. If the integral / is invariant under a [an infinite 
continuous group] Goop depending on arbitrary functions and their 
derivatives up to order <r, then there are p identities among the 
Lagrangian expressions and their derivatives up to order o. Here 
as well the converse is valid. 

Furthermore she has another important result, although it follows easily from 
Theorem II, in our opinion, both because of its importance and the fact that it 
was the key issue motivating her investigation, it could have been set off and called 

Theorem III: 

Given I invariant under the group of translations, then the energy 
relations are improper if and only if I is invariant under an infinite 
group which contains the group of translations as a subgroup. 

Regarding this latter result, she ends her paper with the remarks 

As Hilbert expresses his assertion, the lack of a proper law of 
energy constitutes a characteristic of the “general theory of rela¬ 
tivity.” For that assertion to be literally valid, it is necessary to 
understand the term “general relativity” in a wider sense than is 
usual, and to extend it to the aforementioned groups that depend 
on n arbitrary functions. 27 


The footnote that ends her paper is also of interest: 
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27 This confirms once more the accuracy of Klein’s remark that 
the term “relativity” as it is used in physics should be replaced by 
“invariance with respect to a group.” 

Her result regarding the lack of a proper law of energy applies not just to Ein¬ 
stein’s general relativity theory, but in fact to all geometric theories of gravity: for 
all such theories there is no proper conserved energy-momentum density. 

As a well known textbook expresses it: 

Anyone who looks for a magic formula for “local gravitational 
energy-momentum” is looking for the right answer to the wrong 
question. Unhappily, enormous time and effort were devoted in 
the past to trying to “answer this question” before investigators 
realized the futility of the enterprise. 

Misner, Thorne & Wheeler Gravitatior J® p 467. 


3. The Noether energy-momentum current ambiguity 


Let us begin our technical discussion by first reviewing some background. 

As we will soon show, a well known result is that from a classical field Lagrangian 
density, C^a, d^ipA), via Noether’s first theorem, the translational symmetry of 
Minkowski spacetime leads to a simple formula for the “conserved” canonical energy- 
momentum density 


T» v := 6£C - 


dC 


dd^tpA 


dvVA- 


The divergence of this expression satisfies the identitjQ 


d^ v = -^-d„VA, 
dip A 


( 2 ) 

( 3 ) 


which can easily be directly verified using the definition of the Euler-Lagrange vari¬ 
ational derivative 


SC dC / dC \ 
5ip A ' dip a M \ddp<p A ) ' 


( 4 ) 


The canonical energy-momentum density is a conserved current in the sense that 
“on shell” (i.e., when the Euler-Lagrange field equations are satisfied: SC/5ipA = 0) 
its divergence vanishes. 

The above energy-momentum density has the usual conserved current ambiguity: 


:= T» v + d x U^ x \ 


( 5 ) 


is likewise conserved but defines different energy-momentum values. Essentially, 
one can always adjust by a “curl” a divergence free current. 

'A consequence of assuming that the Lagrangian depends on position only through the field if a- 
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At first thought one might be inclined to follow the rule of sticking with the re¬ 
sults obtained directly from the Lagrangian and the above formula. But sometimes 
the results so obtained are not so suitable physically. 

A simple example in Minkowski space is the Lagrangian density 

C = ~F* V F I1V , F lu/ =d tl A v -d v A lk . (6) 

If one regards this Lagrangian density according to the above paradigm as a function 
of An and d l _ l A u then the above formula leads directly to the conserved expression 

T% = F» a d v A a - -^F^F aP . (7) 

Now the above Lagrangian density (up to a suitable overall scaling coefficient that 
can be set aside for the purposes of this section) can be used to describe Maxwell 
electrodynamics. As is physically appropriate, this Lagrangian density and the field 
equations obtained from it are gauge invariant under the local “gauge” transforma¬ 
tion An —> A^+d^x- However the above canonical energy-momentum density is not 
gauge invariant (nevertheless, as we will see later, it can still be useful physically). 
Naturally, one would generally prefer to have a gauge invariant energy-momentum 
density for electrodynamics. In this particular case there are several ways that 
one can find an alternative to 0: W one can exploit the abovementioned free¬ 
dom m and thereby find “by hand” an “improved” gauge invariant expression, 
(ii) one could regard the Lagrangian as being a function of a one-form A and its 
differential (this is really the proper way to treat electromagnetism—but then one 
needs an extension of the above classical field theory formalism that can accommo¬ 
date form fields; we will discuss such a formalism below), or (iii) one can consider 
that physically any time one has material energy-momentum one must also have 
gravity: the gravitational equations will include an unambiguous gauge invariant 
energy-momentum density])] From a specific gravity theory one gets then a specific 
formula for the energy-momentum density. In this way the ambiguity in the canon¬ 
ical energy-momentum density for any classical field can be entirely removed when 
we consider their gravitational effects. Specifically, in the case of GR, knowing the 
curvature gives, via the Einstein tensor, the symmetric Hilbert energy-momentum 
density. In particular for the electromagnetism example this is 

T^ = F^F va -^F a ^F a p, (8) 

which is no doubt a good choice for the energy-momentum density for Maxwell 
electrodynamics; in fact it is, as we shall see, the same as the energy-momentum 
density that one obtains by regarding the vector potential as a one-form. 

J The existence of a gravitational field will reveal the location of a source with energy-momentum 
even if it has not otherwise been detected. An important example of this is—assuming gravity 
is well described by GR—from astronomical observations it seems that there is a large amount 
“dark matter” in the universe. Clearly, the issue of the proper description of energy for gravitating 
systems can have major consequences for our conception of the physical world. 
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While for electromagnetism one has another criteria (gauge invariance) that 
can be used to arrive at a physically suitable energy-momentum density, for most 
other sources one can only turn to gravity to identify a unique energy-momentum 
density. Hence gravity uniquely detects the energy-momentum density of its sources. 
It may thus seem somewhat ironic that for gravity itself there is no proper energy- 
momentum density. 


4. Pseudotensors 

The Einstein Lagrangian for GR differs from the Hilbert scalar curvature Lagrangian 
by a certain total divergence which removes the 2nd derivatives of the metrico 

2nC E (g afi ,d^g a p) := 

= ,J=gR-d ll (y/=ggf > 'T a f h 5{Z) 

—. 2kG E {gap , dg,gtx0i got0 ) . (9) 

They give the same field equations. The Einstein pseudotensor can be obtained 
from £e using the aforementioned formula for the canonical energy-momentum 
tensor ©: 

:= E - ^-^-d v9aP . (10) 

(Here, following tradition, gothic letters indicate densities.) We have from ([3]) using 
the Einstein equation 

dg(t E u) = jr^-d„g a t 3 = -{2k)- 1 d v g a p = -]■% a0 d v g a p. ( 11 ) 

og a p 2 

Hence, using the vanishing covariant divergence of the material energy-momentum, 
0 = V M (T^) = ^(1^) - = d^ v ) - ^ a0 d v g a p, (12) 

we obtain 

d^„+t£ v )=0, (13) 

a vanishing ordinary divergence, i.e., a conserved total energy-momentum “current”. 
Here we assumed the vanishing of the covariant divergence of the material energy- 
momentum tensor and used Einstein’s equations to obtain an ordinary divergence 
conserved current. But one can argue the other way around, as Einstein did in 
1916. 


k Herer“ i g 7 := \g a> '{dpg\^+d 1 g\p—d\gp 1 ) is the well-known Christoffel/Levi-Civita connection, 

Sff := and k := 87rG/c 4 . 

pc T IP a]’ ' 
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4.1. Einstein, Klein, and superpotentials 

Einstein^ obtained results of the form 


^er^ + t^) 

= 0, 

(14) 


= dxs^u. 

(15) 


= o, 

(16) 


' 9 dd xg ™- 

(17) 


Klein regarded the first three relations as mathematical identities and argued that 
energy-momentum conservation in GR was fundamentally different from that in 
classical mechanics ! 13 l 72 l Einstein did not agree with either of these statements; he 
regarded only m as an identity, which he obtained using a general coordinate 
invariance argument. Now taking the divergence of m using m gives (El), 
which—reversing the computation in the previous subsection—leads to (1121) . In 
this way Einstein showed that the local coordinate invariance identity plus his field 
equations—which are equivalent to m — gives the conservation of material energy 
momentum, without ever having to use any matter field equations. This type of 
argument is referred to as automatic conservation of the source (see MTWp^ section 
17.1); effectively it uses a Noether second theorem type of argument to obtain 
current conservation. Weyl used the same type of argument for the conservation 
of the electromagnetic current in his seminal gauge theory papers p 47 ) 148 ^ whereas 
modern field theory books generally use Noether’s first theorem in connection with 

current conservation.^ 

The identity (flfil) is equivalent to the contracted Bianchi identity, V^G^,, = 0, 
where G^ v is the Einstein curvature tensor. In those days the Bianchi identity , 

i? a/ 3 ^ 7 ] = 0, was not generally well known (it was first used in GR by Levi-Civita 
in 1917H2)). For any Lagrangian constructed out of the metric and its derivatives, 
it is now well known that local diffeomorphism invariance (with 8g^ u = £^g^ v = 
V M £„ + V^) of the associated action leads to a divergence identity: 

[ -j^—dg^d^x = 0 => ^vT~~ — 0- ( 18 ) 

J ogog^ 

In general such identities can involve higher derivatives of the curvature, however 
for the Hilbert scalar curvature Lagrangian of GR this Noether second theorem type 
argument yields a divergence which coincides with the contraction of the Bianchi 
identity. 

By the way, Einstein had been using essentially the set of energy-momentum con¬ 
servation relations El) El) for some years in connection with the (non-covariant) 
“Entwurf” equations that he worked out with Marcel GrossmannPSl However, his 
Lagrangian for that scheme (see CPAEp^ID y 0 i Doc. 2) was not—up to an ex¬ 
act differential—diffeomorphically invariant, so that (HI was in that case a relation 
that selected a preferred set of coordinates. 
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Since for the Einstein Lagrangian © is an identity, one might wonder why 
did Einstein (and others) favor an invariance argument rather than just directly 
calculating it? It turns out that the invariance argument is considerably easier. Let 
us look into this a little further. The detailed form of m can be written out with 
the help of some formulas given by Tolman. 1 ^ From his eqs (87.5) and (89.9) we 
have 

s ,J \ = -Q a ^W\ a + ^g a0 W x aP , (19) 

= -r% + r> CT . (20) 

To show directly that this expression satisfies (flljl) is not short or simple. The only 
published calculation that we know of is in Mpllcr P some 40 years later. 

If were antisymmetric in its upper indices the identity would be trivially 
satisfied. Stated another way, if (TTH) is to be satisfied then there should exist a 
superpotential iF lA „ = U^ A 1„ such that 

‘Z» v + \r v = d x W\. ( 21 ) 

Einstein’s 5 /iA !/ is not antisymmetric; it is not the right kind of potential. A suit¬ 
able superpotential for the Einstein pseudotensor was found over 20 years later by 
Freud 

( 22 ) 

To obtain this he did not follow Einstein’s path. He started with the basics, the Ein¬ 
stein equations, and rearranged them using some formulas from Weyl’s bootf^ and 
some complicated identities he found in Pauli’s 1921 encyclopedia articleP^Si Later 
we will give a simple derivation of Freud’s superpotential using a better technique. 

4.2. Other GR pseudotensors 

The presence of a non-vanishing energy-momentum density necessarily produces 
gravity (i.e., the curvature of spacetime). In curved spacetime the total source 
energy-momentum tensor satisfies m Without the second term we would have 
an expression suitable for integrating to obtain a conservation law. The second term 
represents a local interaction exchanging energy-momentum between the source and 
the gravitational field. To have a good conservation law we would like to rewrite m 
in the form of (11411 for some suitable gravitational energy-momentum density t/V. In 
fact this can be done in an infinite number of ways, and , in all cases, the quantity 
t/V is not a tensor. (For some good overviews of such pseudotensors and their 
properties see Refs. 1451(15118011 381 .1 

Here is a construction. Select an object (referred to as a superpotential) with 
suitable symmetries: il" A At = il^ A l . Now use it to split the Einstein tensor, defining 
a gravitational energy-momentum pseudotensor according to 


2rWA := -2V=?G M „ + d x il' lX v . 


( 23 ) 
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Then Einstein’s equation, G M „ = takes a form (analogous to Maxwell’s equa¬ 

tion) with a total effective energy-momentum pseudotensor as its source: 

dxlPK =2«(n+T^). (24) 

The essential feature of such source expressions is that they are equated to a deriva¬ 
tive of a superpotential in such a way that their divergence automatically vanishes. 
Thus, as a consequence of the symmetry of iff lA I , we have (similar to Maxwell’s 
theory) “automatic conservation of the source”: ^(T'T + V l v ) = 0. This expres¬ 
sion can be integrated to define the total conserved energy-momentum within any 
volume V: 

P v {V) ■■= J (W* + t'*„)d% = J v d x W x v d 3 Z» = ±- j gv W\ l -d>S» x . (25) 

From the volume integral on the left one would expect that the results would be 
highly ambiguous—depending on the choice of reference frame throughout the vol¬ 
ume of interest. However, from the last surface form, one can see that the situation 
is not quite so bad. The result does not depend on the choice of reference frame 
within the volume, it is quasi-local , i.e., it depends on the fields and choice of refer¬ 
ence frame only on the boundary. It should be noted that (for any given a reference 
frame on the boundary) the value of P„(V) is well defined by the above integral. Its 
value, however, comes from a mixture of physics and a quasi-local reference frame; 
still it can be useful if one is mindful of its nature. 

There are some variations on the above formulation. The classical pseudotenso- 
rial total energy-momentum density complexes, T 71 ,, := IT,, + t'T- all follow from 
suitable superpotentials according to one of the patterns 

2nT li v = dxU! lX u, 2 kT i * v = dxUP", 2 kT^ = d aP S) a ^ v , (26) 

where the superpotentials have certain symmetries which automatically guarantee 
conservation: specifically iff lA i; = il[ /,A l v . iP lAl/ = while has the alge¬ 

braic symmetries of the Riemann tensor (this latter form yields a symmetric pseu¬ 
dotensor and, hence, a simpler conservation of angular momentum description, see 
MTWP §20.3). We have already considered the Einstein total energy-momentum 
density which follows from the Freud superpotential m- For completeness we list 
the other well-known ones. The Bergmann-ThompsonP Landau-LifshitzpS Papa- 
petroup^l Weinber^i^ (also used in MTM®) and Mpllei^] total energy-momentum 
complex expressions can be obtained respectively from 


zx BT 

■= g uS Up s , 

(27) 


■= IsI^bt", equivalently Sj^f v := | g\S%* a g af3 g™, 

(28) 


■■=&f b g ab (\g\h mn ), 

(29) 

<~a)iPv 

JJ W 

■■= ^aC\~9\ h ~9 ab {-9 mc g nd + \-9 mn 9 cd )9c d , 

(30) 

U^ A 
• U M v 

:= -Isl V'TV^ = M - dsgp*). 

(31) 
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Here g ab is the Minkowski metric, all indices in these expressions refer to spacetime 
and range from 0 to 3, otherwise our conventions follow MTWP 

People have often looked askance at such pseudotensors, e.g., the above quote 
from MTW and Schrodinger refers to them as “sham” LliS] a s we uojg^ there are no 
doubt two unsatisfactory aspects: (i) which of the many possible expressions should 
one use? (ii) and which quasi-local (in view of (E5l) l reference frame should be used. 
On the other hand one should also be mindful that (a) they do provide a description 
of energy-momentum conservation, (b) they (like connection coefficients) really are 
geometric objects, with well defined values in each reference frame (this issue has 
been rigourously addressed using fiber bundle formulation d 43 | 105 ll201) . 

All of these pseudotensors (except for Mpllcr’s) give the expected total energy- 
momentum values at spatial infinity. On the other hand, none of them give the 
desired positivity of energy for small vacuum regions to lowest non-vanishing or- 
der 11221 however a set of new pseudotensors depending on several parameters with 
this desirable property has been constructedP^ How can one understand the phys¬ 
ical significance of these various pseudotensors? We have found a way using the 
Hamiltonian approach. 


4.3. Pseudotensors and the Hamiltonian 


To see how one can be led to the Hamiltonian one need merely redo the calculation 
of (1^51) as an identity (“off shell”). For some fixed reference frame, with a (constant 
in the present reference frame) vector field Z M inserted we fincfl 


- Z IJ 'P fl (V) ~ - / Z^%V=gd 3 £„ 


IV 


lv l 


Z»yf=-g ( -G\ - T\ ) - —d x (ZW\) 


d 3 ^ 


= f Z^H GR +<£ B GR (Z) = H(Z,V). 
Jv JS=dV 


(32) 


Here Ti GR can be recognized as the covariant expression which, when expressed in 
terms of the appropriate canonical variables, is just the ADM Hamiltonian density 
(i.e., the superhamiltonian and supermomentum), see, e.g., Refs [1160] and MTW^ 
Chapter 21. The expression includes a Hamiltonian boundary term, a 2-surface inte¬ 
gral of B GR {Z) = -Zh\l/2K)'d vX fi {l/2)d 2 S v X: i-e., it is entirely determined by the 
superpotential. The value of the Hamiltonian on a solution is entirely determined 
by this boundary term; the initial value constraints ensure that the Hamiltonian 
density in the spatial volume integral vanishes “on shell” (i.e., when the field equa¬ 
tions are satisfied). In a similar way the value given by any pseudotensor can be 
regarded as the value of the Hamiltonian with a certain boundary termP^ From 
the Hamiltonian variation, as we will discuss below, one gets important information 

'The sign in this expression is dictated by the condition for positive energy determined by the 
Hamiltonian using our local Minkowski signature convention: = ( ~E/c,p). 
















July 28, 2015 0:26 


World Scientific Review Volume - 9.75in x 6.5in 


100GR'Ch4'150726vl page 15 


Gravitational energy for GR and Poincare gauge theories:a covariant Hamiltonian approach 15 

that tames the ambiguity in the boundary term—namely boundary conditions— 
and thereby determines the physical significance of the various quasi-local values. 
The energy-momentum values obtained for the various pseudotensors can all be 
regarded as values of the Hamiltonian with different boundary conditions. 

5. The quasi-local view 

The modern idea, due to Penrose!® 1 in 1982, is that energy-momentum is quasi¬ 
local: i.e., it is associated with a closed 2-surface (while the pseudotensor energy- 
momentum complexes always had this property, its essential importance became 
much more appreciated after this work of Penrose which introduced this convenient 
term). There is a comprehensive review of this topic: Szabados (2009) E3l] The 
many recent works cited in this review show that this is still a topic of considerable 
interest. 

In a brief summary one can find the statement: 

“... contrary to the high expectations of the 1980s, finding an 
appropriate quasi-local notion of energy-momentum has proven to 
be surprisingly difficult. Nowadays, the state of the art is typi¬ 
cally postmodern: although there are several promising and useful 
suggestions, we not only have no ultimate, generally accepted ex¬ 
pression for the energy-momentum and especially for the angular 
momentum, but there is not even a consensus in the relativity 
community on general questions ... or on the list of the criteria of 
reasonableness of such expressions.” 

However if one takes a more specific approach, one can come to a more satisfac¬ 
tory conclusion. In particular the Hamiltonian view quite changes the prospects, 
especially when used along with a gauge perspective. 

6. Currents as generators 

Noether’s work was entirely Lagrangian based. Her results can be taken a further 
step when they are combined with the Hamiltonian formulation. As we will see, the 
Hamiltonian formulation offers a handle on the Noether current ambiguity. 

One key feature can be seen already in Hamiltonian mechanics. A quantity Q 
conserved under the time evolution generated by a Hamiltonian H = H(q,p) is 
more than a just a conserved quantity, it is also the canonical generator of a one 
parameter transformation on phase space (g(A),p(A)) which is a symmetry of the 
Hamiltonian. 



( 33 ) 
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Iii Hamiltonian field theory, the conserved currents are the generators of the as¬ 
sociated symmetry. In particular, the generator of a local spacetime “translation” 
(an infinitesimal diffeomorphism) is the Hamiltonian; energy-momentum is the as¬ 
sociated conserved quantity. Conversely, for spacetime translations, the associated 
Noether conserved current expression (i.e., the energy-momentum density) is the 
Hamiltonian density—the canonical generator of spacetime displacements. As we 
will see, because it can be varied this translation generator gives a handle on the 
associate conserved current ambiguity. The Lagrangian formulation affords no such 
handle, because in terms of Lagrangian variables the translation current is not a 
generator that can be varied. 


7. Gauge and geometry 


For the early history of gauge theory see O’RaifeartaighP^l Briefly, the milestone 
works are Hermann Weyl’s treatments of electromagnetism: Weyl (1918) Weyl 
(1929) pAHl then the generalization to non-Abelian groups by Yang & Mills (1954)Pi] 
and Utiyama (1956, 1959)4^332] Explicitly treating gravity as a gauge theory was 
pioneered by Utiyama 1 139 1 140 1 using the Lorentz group and Riemannian geometry. 
SciamgPd®] also used the Lorentz group but with Riemann-Cartan geometry (i.e., 
non-vanishing torsion). Kibble® put things in their proper place, he gauged the 
Poincare group (i.e., the inhomogeneous Lorentz group, including translations). 

For accounts of gravity as a spacetime symmetry gauge theory, see Hehl and 
coworkers j 46 l 53 ^ ^ 7 Mielkc 81 ^and BlagojevicP A comprehensive reader with sum¬ 
maries, discussions, and many reprints has recently appeared: Blagojevic & HehlP 
For the observational constraints on torsion see Ni (2010).® 

To us it is rather surprising that the idea of regarding gravity as a gauge the¬ 
ory is not better known. Examined more closely, one finds that gravity played an 
important role in the argument used in both of the above mentioned seminal works 
of Weyl, and thus in all of the above—except for the Yang-Mills paper. Further¬ 
more, later in 1974 Yang himself published a papei^® where he proposed a certain 
treatment of gravity as a gauge theory^ 

According to our understanding, properly speaking, GR can be understood as 
the original gauge theory. After all, it was the first physical theory where local 
gauge freedom (in the guise of general coordinate invariance) played a key role[)] 
The conserved quantities, energy-momentum and angular momentum/center-of- 
mass momentum are associated with with the geometric symmetry of Minkowski 


m The aforementioned reader includes a chapter with a critical discussion of Yang’s gauge theory 
of gravity. Recently Yang was asked about his 1974 paper, he said: “I do not believe that paper 
is correct.’t^^ 

n It is true that the electrodynamics potentials along with their gauge freedom were known lon g 
before GR (in fact a Lagrangian which is locally gauge invariant had already been presentecf^^) , 
but this gauge invariance was not seen as having any important role in connection with the nature 
of the interaction, the conservation of current, or a differential identity—until the seminal work of 
Weyl, which post-dated (and was inspired by) GR. 
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spacetime, the spacetime translations and Lorentz rotations, i.e., the Poincare 
group. Furthermore this group is used to classify physical particles according to 
mass and spin. So a local Poincare gauge theory is quite appropriate both geomet¬ 
rically and physically. 

To give a good account, one should also be mindful of the parallel development 
of the closely related concept of a connection in differential geometry. Here we just 
briefly mention that the main ideas were due to Hessenberg, Levi-Civita, Schouten, 
Weyl, Cartan, Ehresmann, and Koszul; for discussions of connections see Nomizuf^ 
Kobayashi & Nomizu -2 and SpivakfA^U 

As we will see in more detail, Riemann-Cartan (with a metric and a metric 
compatible connection, having both curvature and torsion) is the most appropriate 
geometry for a dynamic spacetime geometry theory: its local symmetries are just 
those of the local Poincare group. So in this presentation we will be considering 
the Pincare gauge theories of gravity (PG); GR is included in this class as a special 
case. 

8. Dynamical spacetime geometry and the Hamiltonian 

We will consider geometric gravity theories with both a metric and an a priori metric 
compatible connection. Both curvature and torsion are allowed. The variational 
principles are developed. The Noether symmetries and the associated conserved 
quantities and differential identities are discussed. From a first order Lagrangian 
formalism using differential forms, we construct a spacetime covariant Hamiltonian 
formalism. The Hamiltonian boundary term gives appropriate expressions for the 
quasi-local quantities, energy-momentum, angular momentum and center-of-mass 
momentum, as well as quasi-local energy flux. The formalism easily specializes to 
teleparallel theory and Einstein’s GR. 

The Hamiltonian approach reveals certain aspects of a theory, including the 
constraints, gauges, and degrees-of-freedom, as well as expressions for energy- 
momentum and angular momentum. However the usual ADM approach achieves 
this at a heavy cost: the loss of manifest 4D-covariance. Our alternative approach 
is complementary: a major benefit is manifestly 4D-covariant expressions for the 
quasi-local quantities: energy-momentum and angular momentum/center-of-mass 
momentum. 

8.1. The main ideas 

The Hamiltonian for physical systems and dynamic spacetime geometry generates 
the evolution of a spatial region along a vector field. It includes a boundary term 
which determines the boundary conditions and supplies the value of the Hamilto¬ 
nian. The Hamiltonian value gives the quasi-local quantities: energy-momentum 
and angular-momentum/center-of-mass momentum. A spacetime gauge theory 
perspective identifies suitable geometric variables. We found a certain preferred 
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Hamiltonian boundary term. The Hamiltonian boundary term depends not only 
on the dynamical variables but also on their reference values; they determine the 
ground state—the state with vanishing quasi-local quantities. To determine the 
“best matched” reference metric and connection values for our preferred boundary 
term we propose on the boundary 2-surface: (i) 4D isometric matching, and (ii) 
extremizing the energy. 

8.2. Some comments 

Before we begin our technical discussion of our work, let us make a few general 
comments. We work in 4D spacetime, but most of this can be extended to other 
dimensions in a straightforward fashion (except for the reference construction). The 
class of dynamical Lagrangians we will consider does not allow for any derivatives 
of curvature or torsion. Our concerns are entirely classical. 

We focus here on Riemann-Cartan spacetime geometry (i.e., spacetimes with 
a metric and a metric compatible connection, having both curvature and torsion) 
and the PG; our general analysis can be specialized both to Riemannian geometry 
(vanishing torsion) and teleparallel geometry (vanishing curvature); it includes GR 
and the teleparallcl equivalent of GR as two special cases. Here we assume a metric 
compatible connection; elsewhere we will present the generalization which includes 
non-metricity. The extension to non-metricity and the special case of teleparallcl 
geometry each offer further insight into gravitational energy; we believe those in¬ 
sights are best appreciated when compared to the results presented here for the 
Riemann-Cartan geometry with the PG. 

9. Differential forms 

In this work we mainly use differential forms )Hl 56 l 142 l The reader may wonder why 
we use this less widely familiar idiom. The simple brief explanation is that they 
have some qualities that are technically very convenient for our needs. Differential 
forms are multiplied using the (graded Grassmann) wedge product. They can be 
differentiated using d, the exterior differential , a graded derivation which enjoys the 
property d 2 = 0, so a differential equation da = /3 has the integrability condition 
0 = dj3, furthermore d{3 = 0 => /3 = da, at least locally. The integrals of forms 
satisfy the general boundary theorenf] 



(34) 


Also, they are well suited to representing interacting physical fields, especially gauge 
fields, and, as we shall see, they give a succinct representation of the main geometric 
objects: connection, curvature, coframe, torsion, etc. Moreover, as will be explained 

°This generalization of the fundamental theorem of calculus is often referred to as the generalized 
Stokes theorem. Special cases include the Ostrogradsky-Gauss and Stokes theorem of vector 
analysis. 
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below, they are quite convenient for the essential 3+1 spacetime decomposition of 
derivatives that is needed for a dynamical Hamiltonian formulation.^ 

Regarding notation, here the contraction (or interior product) with a vector 
field is denoted by ix<x(-, • ,•■•) : = a (X, (some authors use the notation of 
left contraction: Xj a). The Lie derivative on forms is given by £x = i.xd + dix , 
it has the nice property d£x = £xd. We are concerned here with the case of 
4-dimensional spacetime, which has a local Minkowski structure, having a metric 
with Lorentz signature. The metric determines the unit volume 4-form 77 with 
components 770123 = \/\g\ which is used to construct the Hodge 

dual that maps fc-forms to (4 — fc)-forms. 

From the coframe one can construct a basis for fc-forms $ a P--- := /\d^ /\. .. 

and a useful dual basis 77 “ /3 -" := They are related by various identities, 

especially 

d p A = Sfav + S^rjvx + <5£?7a m , (35) 

d af3 A r) X + 77 ^ + S 77 ^, (36) 

d p A 77 ^ = 5 p ^ -5 p r] v . (37) 

Maxwell’s electrodynamics is a good example of the utility of differential forms. 
Charge identifies the charge-current 3-form (density): 

<200 = f J- (38) 

Jv 

Charge is conserved: 

dJ = 0 =+ J = dH. (39) 

The electromagnetic field is represented by a 2-form F. An elementary way to 
see why this is appropriate is to examine the motion of a point test charge. One 
should begin with kinematics in Minkowski space. Consider the motion of a point 
particle as a function of proper time: x p = a+(r). The 4-velocity v p := dx p /dr has 
constant magnitude: v p v^ = —c 2 so the 4-acceleration is Lorentz orthogonal to the 
4-velocity. Hence the 4-force must be orthogonal to the 4-velocity. Consequently 
the 4-force must depend on the velocity. The simplest case is for a 4-force linear in 
the 4-velocity. Thus the simplest dynamical law has the form 

= qF^v v , (40) 

where p M := mv^ is the 4-momentum, q is a coupling constant and F is some 
tensor field which is antisymmetric , i.e, it is a 2-form. The Lorentz force law of 
electrodynamics has this form. The force law identifies a certain field strength 2- 
form F which includes the electric and magnetic fields. Conservation of magnetic 
flux through a closed 2 surface S = dV gives 

0 = / F = I dF, => dF = 0 +=+ F = dA , (41) 

Js Jv 
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and there are local gauge transformations: A —► A + d\- The vacuum constitutive 
relation is H = *F/Zq (Zq is the vacuum impedance). This covariant formulation 
for Maxwell’s electrodynamics is valid for all dynamic geometry gravity theories and 
does not depend upon using a particular set of units, for a detailed, comprehensive 
and instructive presentation see Hehl & ObukhovP^ 


10. Variational principle for form fields 


Why do we use variational principles!^ The answer is pragmatics: because they 
work. With appropriate symmetries they give consistent interacting field equations 
along with conserved Noether currents for all the desired quantities. As far as 
we know, all the known good dynamical evolution equations for the fundamental 
interacting classical field theories have a variational formulation. 

In the usual formulations most dynamical fields satisfy second order equations. 
We refer to such formulations as second order. 

Let ip A be some kind of vector field. The label “A” stands for some collection of 
indices, e.g., spinor, spacetime, isospin. Allow ip A to also be a differential form of 
rank / where / = 0,1, 2, or 3, e.g., ip A = \<p A u Ad v = ^ipfjdx 1 A dx J for / = 2. 

The Lagrangian density is a 4-form: 

£ = C{p A ,dp A ). (42) 


Note that there is no explicit appearance of the coordinates x l or the coordinate 
partials dp, dp A is an (/ + l)-form which geometrically includes partial derivatives 
of the components of ip A , but only in an antisymmetric combination. (Here we 
explicitly consider just one /-form field. The generalization to include several fields 
of different grades, is straightforward.) 

Our convention is to vary fields off to the left (other conventions would differ 
only by some signs). The variation of £ is thus 

SC = A 8^* + Sv " A W' <43) 

This implicitly defines dC/ddip A as a (3 — /)-form and d£/d<p A as a (4 — /)-form. 
Next, interchange the order (i.e., Sd = dS) to get 


S£ = dS<p A A 


= d 


d£ 
dd<p A 
d£ 


Sip A A 


d£ 

dip A 


5vA A aS^r) “ ( - 1) '^ 4 A d (aS^O + V 


A 


dC 

d<p A 


(44) 


(Upon integration over some spacetime region this last step is “integration by parts”, 
with the total differential term becoming a boundary term.) From the above it 
follows that the basic variational relation 

SC = d(Sip A A pa) + 5ip A A 

Oi p A 


(45) 
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can be regarded as implicitly defining the conjugate field momentum and the Euler- 
Lagrange variational derivative , which have the respective explicit definitions 


dC 

PA ddp A ’ 


(46) 


SC dC ( dC \ 

5p A dp A ^ \ddp A ) 

A small price for using form fields is the appearance of occasional sign factors like 

C := (-1)/. 


10.1. Hamilton’s principle 

Our first application of (1451) is Hamilton’s principle {the principle of least action). 
Let the action within a region U be given by S := f v C. Then 

SS= [ SC = [ d{Sp A Ap A ) + 6p A = [ Sp A A y^t+S p A Ap A . (48) 
Ju Ju Sp A Ju Sp A J du 

Now require the action S to be extreme (i.e., SS = 0) for 6<p A vanishing on the 
boundary of U. This yields the field equation SC/Sp A = 0. 


10.2. Compact representation 


For a compact general discussion it is convenient to suppress the field component 
index.® (This could be represented in matrix notation; our basic fields ip and their 
differential could be regarded as row vectors.) The Lagrangian then has the form 
C = £(p,d<p) and the variational scheme proceeds as 


dC 

SC = dSp A ——I- Sip A 
odp 


dC 

ddp 


= d ( Sip A ^ + Sip A 


, / dC \ dC 
-qd I - 1 -l- 


\ ddip J dp 

hence the key Lagrangian variational identity takes the form 


(49) 


SC = d(Sp A p) + Sp A —. (50) 

dp 

In this succinct alternative to (1451) the conjugate momentum and the variational 
derivative can be regarded as form-valued column vector fields. 


11. Some simple examples of the Noether theorems 

Here we present simple examples of Noether’s two theorems.® Later we shall use 
the same types of arguments in more complicated situations. 
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11.1. Noether’s first theorem: energy-momentum 

Further applications of the basic variational identity (H5l) or (15U1) yield the Noether 
theorems. Their applications to physical systems are discussed in many works, e.g., 


Konopleva and Popov.^ Here we introduce our particular use of them using two 


specific important cases. 

Noether’s first theorem states that for a constant parameter symmetry there is 
a conserved current. 

For our concerns the most important example is the conservation of energy- 
momentum. As our specific relevant simple case exemplifying the argument, we 
specialize in this subsection to Minkowski spacetime, which is homogeneous and 
thus naturally has a geometric symmetry under translations. Dynamically let us 
assume symmetry also of the action under constant translations. The symmetry 
depends on a continuous parameter; it is sufficient to consider the infinitesimal 
case. Geometrically an infinitesimal translation corresponds to a constant vector 
field Z. Under such a transformation the change in the components of form fields 
is given by the Lie derivative. We have: 


Ap = -£ z p = -(diz + iz d)p, 
AC = —£zC = —diz£- 


(51) 

(52) 


Equation (15U1) under these specific variations (i.e., replacing the general S by these 
specific changes) should be an identity (since the Lagrangian C depends on the 
position only through the fields <p). Rearranging leads to 



(53) 


Because of (1521) the l.h.s. of eq. (1531) is a total differential—a total differential 
which moreover vanishes if the Euler-Lagrange field equations arc imposed: 



(54) 


This identifies a conserved current density (a 3-form with vanishing differential on 
shell , i.e., when the field equations are satisfied) called the generalized canonical 
stress energy-momentum density (3-form): 


T(Z) := i z C - £ z p Ap. 


(55) 


For a zero-form field it takes the special shape T a fl Zh rj a where (with L = *C) 



( 56 ) 


which is the well known expression mentioned earlier Q . 
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11.2. Noether’s second theorem: gauge fields 

Note: For the rest of the present section there is no need to restrict the spacetime 
geometry in any way. Our considerations apply quite generally. 

Now we want to consider invariance under local gauge transformations: 

Ap = a p pT p , (57) 

where the a p are position dependent parameters and the T p are the matrix genera¬ 
tors (in our representation the fields are on the left and the matrix on the right) of 
the gauge group. Replace dp by the gauge covariant differential: 

Dtp := d<p + A p A pT p , (58) 

containing a certain compensating field, the gauge vector potential (a.k.a. the gauge 
connection one-form): A p = A p jdx J , which has a special non-homogeneous gauge 
transformation: 

AA P = -Da p := -(< da p + A q C p qr a r ), (59) 

where C p qr are the gauge group structure constants: [T q ,T r ] = C p qr T p . Then 
A (Dip) := (A D)p + DAp 

= -(DoF) A pT v + D(a p pT p ) 

= a p (Dp)T p . (60) 

Thus Dp transforms just like p. 

Rather than starting with a Lagrangian 4-form of the type £ = £(pp, dp, A p , dA p ) 
and then discovering that the variables A p , dp, dA p can only appear in nice covariant 
combinations, let us proceed more covariantly, beginning with the Lagrangian 4- 
form 

£ = £{p,Dp,A p ,F p ), (61) 

where F p is the field strength or gauge curvature 2-form: 

F p := dA p + i C p qr A q A A r . (62) 

Since A p is still allowed to appear independently in (l6ll) . there is no loss of generality. 
The variation of £ (liill) is 

r d£ d£\ S£ 5C 

\ dDp dF p J dp 5A p y ’ 

Now we are set for an example of Noether’s 2nd theorem—that for each local 
invariance there is a differential identity. 

Assume that £ (El) is invariant under the special changes Ap, A A of eqs El 
[59]). From (l63l) we then have the identity 

0 = d (cCyTp A - DcC A AL) + a P v A T p A _ ° a P A AL. (64) 
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dC 


The second term in the parenthesis may be rewritten as —d[a p Appp) + aPD-gp^, 
then using d 2 = 0 gives 


0 = d 


ipT p A 
= Da p (ipT p A 


dC 

dDip 

dC 


dDip 


D^-\ 

dF p J 

dC 

dF p 


a p ipT p A — Da p A SL 

dip 


SA p 


D 


SC 

SA p 


+a p 


D ipTp A 


dC ^dC\ „ SC 
+ D dF^J +<pTpA ~ 


(65) 


( 66 ) 


dDip ' dF p J ' -r-v' 5ip ~ ■ 

For a local symmetry the quantities a p and Da p are pointwise independent; their 
coefficients must vanish separately. The coefficient of a p identifies a Noether I type 
conserved current: 


Jp :— A 


which satisfies the “conservation” law 


dC 

dDip 


D 


dC 
dF p ’ 


SC 


DJ p = -ipTp A —. 


(67) 


( 68 ) 


The r.h.s. vanishes “on shell” (i.e., when the field equations are satisfied). 

From the coefficient of DaP we obtain an algebraic identity relating the Noether 
I current to a variational derivative: 

Jp ~ Hap’ (69) 

thereby the Noether I current conservation becomes a differential identity 

(70> 

between the variational derivatives. Note that to obtain these results there is no 
need for the explicit form of the field equations. 

Another way to argue is to replace the last term in (R)5l) with a total differential 
minus a compensating term, bringing that relation into the form 


0 = d 


ipTp A 


dC 

dDip 


D 


dC 
d Fp 


SC 

5A p 


ipTp A 


5C^ d SC_\ 


Sip 


SA p J ' 


(71) 


If one integrates this over any region the total differential term gives rise to an inte¬ 
gral over the boundary. To have a vanishing value for all possible gauge parameters 
with small support, the coefficient of the gauge parameter everywhere within the 
region and the coefficient of the gauge parameter everywhere on the boundary must 
both vanish identically. This again yields m and 


ipTp A 


dC 


dDip 

which is equivalent to (1681) with (1691) . 


D 


dC SC 


dFP SA p 


= 0 , 


(72) 
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11.3. Field equations with local gauge theory 


It should be noted that the Noether invariance argument yields the differential 
identities just found involving the Euler-Lagrange expressions without any need to 
have the explicit form of the Euler-Lagrange expressions. Of course if one explicitly 
computes the Euler-Lagrange expressions one could go on to verify these identities 
directly. Furthermore, if one has the Euler-Lagrange expressions one could (proba¬ 
bly not so easily) directly discover such identities, even if one was not aware of the 
local symmetry. 

To compute the field equations the explicit variations 


5 Dp = DSp + SA P A pT p , (73) 

SF p = dSA p + C p qr A q A SA r = DSA P , (74) 

are needed. The variation of the Lagrangian 4-form C (IH3l) is 


„„ dC _ dC 

SC = SDp A A ^ 


■ 5F P A 


dC 

dFP 


SA p A 


dC 

8Ap 


= ( D5(p + SA P A pT p ) A 


= D ( 6<p A 


+SA P A D 


dC 

dD^> H 
dC 


dC _ dC dC 

__ + ^ A _ + ^A —d 

+ Sip A + ^ 


SA p A 


dC\ 

SA P A t— 

dFP J 

dC m dC \ 
+ (pTpA dD^J- 


dDp dp 


dFP dAP 


dC 

dAP 


(75) 


Comparing the explicit form of SC/5A p found here with (15511671) shows that (1691) 
means 


dC 

dAP 


= 0. 


(76) 


Thus local gauge invariance means: no explicit dependence on the gauge potential; 
all dependence on A p comes through Dp and F p . Furthermore if one makes the 
usual minimal coupling assumption, 


C = C A (A p ,F p )+C v (p,Dp,A p ), (77) 

the identities (EG!) apply separately to each term. Hence, the Lagrangian 4-form 
must have the simpler form 


C = C A (F p )+C v (p,Dp), 

and the differential identity (1701) becomes the two identities 


D S Ai s 0 , 

SAP 


n SC v _ 5C V 

D-r- = —pJ- v A 
SAp ^ p Sp 


The first relation is explicitly 


0 = £>' 


_ = _ pqrir A A 

dFP - ,p dF'- ■ 


(78) 


(79) 


(80) 
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The latter is a kind of gauge current “conservation”, as the r.h.s. vanishes since 

on shell. In more detail, this gauge current “conservation” relation has the form 


(82) 


0 = DJ P = dJ p - A q C r qp A ,/ r . 


Thus it has some similarities to the vanishing covariant differential of the material 
energy-momentum. Just as in that case, one can rearrange the field equation to 
obtain a conserved gauge pseudocurrent. 

We have gone into considerable detail in this relatively simple example. We have 
done this to prepare the reader, because we are going to use a very similar argumen¬ 
tation in connection with the rather more complicated case involving gravity and 
the dynamic spacetime geometric symmetries associated with energy-momentum 
and angular momentum. It will be seen that almost every step used in our later 
argument and every expression has an analogue with what we have done in this 
subsection. 

12. First order formulation 

In this section we discuss the general formulation of the first order formalism; the 
spacetime geometry has no restrictions. 

We proceed from the action principle. Any action principle can be rewritten 
in an equivalent form, which (following, e.g., ADAfO and Kuchai^Sl) we refer to as 
first order, this is the most convenient form for our purposes. Here we present a 
simple argument (essentially the same Legendre transform idea as is used in classical 
mechanics to construct the Hamiltonian) which is applicable to a large class of 2nd 
order Lagrangians. 

Given a 2nd-order Lagrangian 4-form C(p, dip) we define its associated canonical 
momentum in the usual way: 




(83) 


Next we define a 4-form by 


A (ip, dip , p) := dip A p — C. 


(84) 


Now consider the variation of A: 


SA = 5(dp) Ap +dp A Sp — SC 



= dp A 5p — Sp A ——. 

op 


( 85 ) 
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Which shows that A can be regarded as a function only of vb-PO 

This construction takes one from the usual 2nd order Lagrangian to a first order 
type of variational principle: 

£ lst (vb dip,p) = dip Ap — A(ip,p), (86) 


where ip and p are now regarded as being independent variables and are varied 
independently. Varying (18611 gives 


SC lst 


5C lst 


Sdip Ap + dp A Sp — 6 A 
d(Sip A p) — <;5ip A dp + dip A Sp — Sip A 
<5£ lst <5£ lst 

d(5ip A p) + Sp A —-1--— A Sp. 

dip op 


dA 

dip 


dA 

dp 


A Sp, 


(87) 


(Note: we find it more convenient to vary our momentum fields p off to the right. 
This reduces a little the number of appearances of the sign factor and merely 
amounts to a sign convention on the definition of dA/dp.) 

Using independent p and ip variations gives a pair of 1st order field equations 
for the differentials of the fields: 


5C lst 

Sip 



0 = 


5C lst 

Sp 



( 88 ) 


13. The Hamiltonian and the 3+1 space-time split 

Here we introduce the Hamiltonian and the space-time split. In this introduc¬ 
tory subsection we use for motivation some well-know elementary expressions in 
Minkowski spacetime. In the subsequent subsections the spacetime geometry is quite 
general. 


A key feature of the canonical Hamiltonian formulation is that the field equations 
are decomposed into two sets: the initial value constraint equations and the dynamic 
equations. A familiar example which illustrates many of the ideas is Maxwell’s 
vacuum electrodynamics^ The 4-covariant equations were given earlier ®ED: 
d * F = ZqJ, dF = 0, or in tensor index form 

d,(V^gF^) = ZoV^J", d [a F/ v] = 0. (89) 

p The procedure becomes technically somewhat more complicated if (f83l) cannot be inverted for dip 
in terms of p, p. In that case one must introduce some additional variables that appear in A only 
algebraically and thus function as Lagrange multipliers introducing some algebraic constraints. 
We will not go into such complications in our general development here. Later in our treatment of 
Einstein’s GR we will see a concrete example. Examples of how this has been dealt with in field 
theory in practice can be found in.HnHHHH 

q Zo is the vacuum impedance, £o = (Zqc)~ 1 is the vacuum permittivity, and po = Zqc -1 is the 
vacuum permeability. Here we are taking for simplicity poeo = c -2 = 1 in relativistic space-time 
units. 
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They split (in Minkowski spacetime) into the familiar initial value constraints (spa¬ 
tial projections, with no time derivatives): 

V • E = —, V • B = 0, (90) 

£o 

and the time projections, a pair of dynamic equations: 

B + VxE = 0, VxB = /r 0 J + E, (91) 

which contain the first time derivatives of the dynamical fields linearly. 

The canonical Hamiltonian form of these equations is in terms of the 4-vector 
potential (which satisfies F = dA and splits into the scalar and vector potential). 
The familiar vector form is 

constraint V • E = —, (92) 

£o 

dynamic A = — E — V4>, E = —V x (V x A) + ^oJ- (93) 

The scalar potential field appears here, but it has no evolution equation; it can be 
chosen freely. This gauge freedom affects the evolution of an “unphysical” part of 
the vector potential. Considering this along with the constraint on E (1^1) . one finds 
that the electromagnetic field has two physical degrees of freedom. 

13.1. Canonical Hamiltonian formalism 

The canonical Hamiltonian formalisnPO is of interest because it clearly reveals the 
constraints, gauges, and degrees of freedom, as well as the total energy momentum— 
and it offers a practical way to numerically calculate solutions. 

The dynamical theories of interest all have constraints. The canonical formal¬ 
ism for constrained Hamiltonian systems was developed mainly by Dir at® ^1 anc [ 
BergmarniPEI For a general discussion see, e.g., Hanson, Regge, & Teitelboimpl 
Sundermeyer. 128 >1221 Rosenfeld^^ seems to have been the first to consider a Hamil¬ 
tonian approach to general relativity, but this early work was not followed up. As far 
as we know Pirani, Schild and Skinner^^ were the next to address the issue. Dirac 
gave a rather complete treatment in 1958pM3 The treatment by Arnowitt, Deser 
& Misner (ADM) 1 has come to be regarded as the standard. For a basic discussion 
see MTW^ Ch 21 or Isenberg and NesterP^ For some critical comparison, see 
Kiriushcheva and Kuzmin.Going beyond Einstein’s theory, a remarkable “if con¬ 
straint” formalism was developed for the PG by Blagojevic and NikolicPto deal with 
a conditionally degenerate kinetic Hessian. They use the Dirac type of approach; 
so to construct the Hamiltonian one must first find the primary constraints—which 
depend on the conditional degeneracies of the Legendre transformation. The “if 
constraint” technique is a marvelous way to manage the technicalities involved in 
constructing the Hamiltonian. In our first order approach, in contrast, one can read¬ 
ily formally construct the Hamiltonian and the Hamiltonian equations, however (in 
line with the principle of the “conservation of difficulties”) a suitably adapted ver¬ 
sion of the “if” constraint technique will still be needed when one actually tries to 
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solve the dynamical and constraint equations. The first order approach as used in 
the covariant Hamiltonian formalism allows one to investigate the general formalism 
and, in particular, to find covariant expressions for the “conserved” quantities while 
postponing dealing with such technical details. 

13.2. The differential form of the space-time decomposition 

Note: for the rest of this section, the spacetime geometry is quite general. 

A feature of this standard approach is the loss of manifest 4-covariance. Now a 
Hamiltonian formulation essentially requires that the time derivatives to be singled 
out from the spatial derivatives, so in this sense it cannot be truly 4-covariant. The 
usual approach, however, departs far more from 4-covariance than is necessary, all 
the indices are 3+1 projected, leading to much extra bookkeeping. In the ADM 
approach the spacetime metric is replaced by the spatial metric and the lapse and 
shift. However only the derivatives d M really need to be projected. Since interac¬ 
tion fields are one-form fields, this means decomposing the exterior differential d, 
decomposing d will inevitably involve decomposing the differential form. One of the 
reasons for using differential forms is in how nicely they decompose in this fashion. 

Begin from the basic first-order form Lagrangian ( 1861 ) . Its variation ( 1871 ) identi¬ 
fies the first order Euler-Lagrange expressions ( 1881 ) . According to Hamilton’s prin¬ 
ciple, the first order Euler-Lagrange expressions should vanish. This gives us our 
first order field equations. 

We want to extract the “time derivative” of p and <p, i.e., the change with 
respect to an evolution parameter (which we refer to as time) as seen by observers 
who move along some fixed congruence of worldlines. This change is given by the 
Lie derivative in the direction of the (fixed) vector field Z tangent to the congruence: 
i.e., dt '■= £z- The Lie derivative on the components of differential forms is given 
by a simple neat expression: £z = diz + izd. Using this, from the differential we 
can extract the “time” derivative: 

i z dj3 = £ Z P - di z /3 = /3 - dizff (94) 

(The congruence need not actually be timelike. Indeed, what we are doing here does 
not require a metric tensor. Even when one has a metric whether the vector field is 
“timelike” is not an important issue, our whole Hamiltonian formalism is linear in 
the spacetime displacement vector field Z, so by considering the difference between 
two timelike displacements one could get a spacelike displacement. A metrically 
timelike displacement is important when one actually tries to find a physical solution 
to the equations; for evolution one wants hyperbolic equations.) 

The description of “time” also includes the idea of “instants of time”. Geomet¬ 
rically this is a set (foliation) of (non-intersecting) 3-dimensional hypersurfaces (we 
usually think of them as being spacelike). Locally, we can always choose adapted 
coordinates: a+ = {t,x k }, where k = 1,2,3 so that the spacelike hypersurfaces are 
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£ t with t = constant. With respect to these adapted coordinates Z is the directional 
derivative in the time direction: Z = <9 t . Note that izdt = 1. 

From these considerations we are led to define the “time” and “space” projec¬ 
tions of differential forms. We use the notations 


a := izoi , a:= a — dt A a, 


(95) 


to indicate the “time” component and the “spatial” part of a form. These projec¬ 
tions have simple expressions in terms of adapted coordinates 0 
Thus a general form decomposes according to 


a = dt A a + a. (96) 

In our formalism %z and t are thought of as freely (except that izdt = 1) chosen 
covariant fields, then the decomposition of a is essentially covariant. With this 
notation the differential decomposes according to da = dt A da + da. From HMD 

da = d — dd; (97) 


thus we can extract the part with the time derivative. 

It is convenient to decompose the differential operator d itself. In view of the 
adapted coordinate expression d = dx M A = dt A d t + dx k A dk, we define the 
decomposition as 

d=dtAd + d, with d:=£z■ (98) 


Now we can examine the spacetime decomposition of our first order field equa¬ 
tions (|88l) . We first consider the time projections, which include all the time deriva¬ 
tives: 


S£ lst \ 
Sip J 


( 8 C lst 
V S P 


<£-J7)-(D)=o. 


(99) 

( 100 ) 


they are the dynamic equations for p and ip. In order to use these equations to 
evolve p, we generally need to know p and ip, which normally are provided by 
the initial value constraints: the spatial restriction of (1M1) : 



( 101 ) 

( 102 ) 


If these two equations can be solved for p and ip all is well and good (in that case 
the two equations are, in Bergmann’s terminology, second class constraints). They 
then define p and ip for all time as functions which depend on ip, p, dp and dp. How 

r We have long used this type of decomposition beginning with(£HlMl sefh^J for a similar technique. 
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to proceed for the case where these quantities cannot be found from the constraints 
is best understood from concrete examples. For our purposes of this work, the 
already-discussed Maxwell electrodynamic example is sufficient. In that case there 
is some undetermined gauge freedom. 


13.3. Spacetime decomposition of the variational formalism 


We decomposed the equations. One could decompose the Lagrangian or its vari¬ 
ation. Our approach easily relates these alternatives, as can be seen from the 
following. 


£ lst = dip Ap — A 
SC lst = d(Sp A p) 



^r lst 

+i Sp A 

Sp 

<5£ lst 

+ r 

> 

<^i 

op 


extract 

[ 


SC lst 

sc lst 

Sp 


3+1 


3+1 


3+1 


3+1 


= (ip — dip) A p — (A + gdp A p) 


— (Sip A p) — d(Sip A p + qSp A p) 


+Sp A 


SC 


1st 


Sp 


- (SC lst \ 


-<r 


fSC lst \ fSC lst \ . 


extract 

(!f). (t) 

V sp j V Sp j 

(f?) 

V sp j y Sp j 


(103) 


(104) 


(105) 

(106) 


For our objectives here we will not need to use this projection into the space and 
time parts of form expressions very much. Our intention here was to include enough 
of the details so that the reader can have some confidence that this formalism can 
yield a proper Hamiltonian description. From what we have discussed, it can be 
seen that dynamical equations in this first order covariant form already contain 
both the constraint and dynamical evolution equations. It should be noted that 
the first line of the above set of relations ( 11031 ) shows how the Hamiltonian can be 
simply extracted from the first order Lagrangian. 


14. The Hamiltonian and its boundary term 

In this section we establish some of our main formal results concerning the covariant 
Hamiltonian and its boundary term. The geometry is quite general. The energy, as 
well as the other conserved quantities, of a physical system can be identified with 
the value of the Hamiltonian. In particular for a gravitating system the associated 
Hamiltonian is proportional to the field equations which vanish on-shell. Therefore 
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the corresponding conserved quantities are determined by the Hamiltonian boundary 
term. The choice of Hamiltonian boundary term is associated with the specific 
boundary conditionl 1 


14.1. The translational Noether current 


The action should not depend on the particular way points are labeled. 
Thus it should be invariant under diffeomorphisms, in particular, infinitesimal 
diffeomorphisms—a displacement along some vector field Z. From a gauge theory 
perspective such displacements are a “local translation”. Under a local translation 
quantities change according to the Lie derivative. Hence, for a diffeomorphism in¬ 
variant action the key variational relation (1ST1) should be identically satisfied when 
the variation operator 6 is replaced by the Lie derivative £z {= diz + izd ): 

<i£ lst ri£ lst 

diz£ lbi = £z£ lst = d(£z<p A p) + £zP A —r-1-j— A £ zP- (107) 

dip op 

This simply means that £ lst is a 4-form which depends on position only through 
the fields tp, p. (According to our understanding this is only possible if the set of 
fields in £ lst includes some dynamic spacetime geometric variables: gravity.) 

From (11071) it directly follows that the 3-form 

H(Z) := £ z pf\p-i z C lst (108) 


satisfies the identity 


<5£ lst <5£ lst 

— dH(Z) = £ z p A—-1--—A £ zP\ (109) 

op op 

thus it is a conserved “current” on shell (i.e., when the field equations are satisfied). 
Substituting (1M1) into (11081) gives the explicit expression 


'H(Z) = d(izp A p) + <,iz£ A dp + <;dp A izP + iz A, (HO) 


thus this conserved Noether translation current can be written as a 3-form linear in 
the displacement vector plus a total differential: 

H(Z) =: Z^H/j. + dB(Z). (Ill) 

Compare the differential of this expression, dH(Z) = dZ^ KH^ + ZhdH^, with (11091) : 
equating the dZ ^ coefficient on both sides reveals that 

SC lBt 5C lst 

Z^H^ = -i z p A ■ +<; Sp A izP' ( 112 ) 

Thus, as a consequence of local diffeomorphism invariance, vanishes on shell ; 
hence conservation of the translational Noether current (11091) reduces to a differ¬ 
ential identity between Euler-Lagrange expressions. This an instance of Noether’s 
second theorem, and, moreover, it is exactly the sort of case to which Hilbert’s 
remark regarding the lack of a proper energy law applies. 
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From the above it follows quite generally, as we remarked earlier in the special 
case of GR ( 1551 ) . the value of the conserved quantity, —P(Z,V), associated with 
a 3D region V is determined by a 2-surface integral over the boundary, i.e., it is 
quasi-local: 



(113) 


For any choice of Z this expression defines a conserved quasi-local quantity. What 
do these values mean? As we shall see in detail later, for a suitable time- 
like(spacelikc) quasi-translation displacement on the boundary the expression de¬ 
fines a quasi-local energy (momentum), and for a suitable quasi-rotation(boost) it 
defines a quasi-local angular momentum(center-of-mass momentum). However it 
must be noted that, like all other conserved currents, the translational current is 
likewise subject to the usual ambiguity: one can add by hand the differential of any 
2-form and still have a conserved current. But that amounts to being able to adjust 
B freely, consequently one could obtain almost any quasi-local value. The Hamil¬ 
tonian perspective brings this freedom under physical control. As we shall show, 
the first order translational current 3-form is something more: it is the generator 
of local diffeomorphisms , i.e., the Hamiltonian. 

14.2. The Hamiltonian formulation 

From the first order field equations ( 1551 ) . by contraction with a “time evolution 
vector held” Z we get a pair of Hamiltonian-like evolution equations for the “time 
derivatives”: £z l P, £zP- A key identity involving these time derivatives is revealed 
by comparing two relations. Consider the projection of the Lagrangian 4-form 
iz£ lst , which from (110811 is just £zp Ap — H{Z)] its variation is 


Siz£ lst = S{£z<p Ap)~ S'H(Z) 


= 8(£zv) A p+ £zp A Sp — S'H(Z) 
= £z5p Ap + £zp A 8p — 8H(Z ) 


= £z(5ip A p) — Sip A £zp + £zp A 5p — 5H{Z). (114) 

Compare this with the projection of <5£ lst ( 1571 ) along Z: 



Since Z is not varied the two relations are identical: 5iz£ lst = izSC lst ; conse¬ 
quently, 
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The last term vanishes “on shell”. This relation identifies the Noether translational 
current 'H(Z) as the Hamiltonian 3-form (i.e., density), as the following considera¬ 
tions show. The integral of 'H(Z) over a 3-dimensional region, 


H(Z,X) := f H(Z), (117) 

J s 

is the Hamiltonian which displaces this region along Z, since the integral of its 
variation: 


5H(Z, E) = J SH{Z), (118) 

yields, from (11161) . “on shell” (then the last bracketed term vanishes) the Hamilton 
equations: 


£zP 


SH(Z , E) 
Sp 


£zp = 


5H(Z, E) 

Sip 


(119) 


if the boundary term in the variation of the Hamiltonian vanishes. In this case 
that means when Sip vanishes on <9E. Technically, the variational derivatives of the 
Hamiltonian H(Z,T,) displayed in (11191) are only defined for variations satisfying 
this boundary condition. In other words this Hamiltonian is “well-defined”, i.e., 
functionally differentiable , only on the phase space of fields satisfying the particular 
boundary condition Sip\gs = 0. 


14.3. Boundary terms: the boundary condition and reference 

In some important cases the fields of physical interest do not satisfy the boundary 
condition naturally inherited from the Lagrangianp^l this happens in particular for 
the spacetime metric of an asymptotically flat region. A modified formulation is 
needed to deal with this. 

One alternative is to modify the Lagrangian 4-form itself by a total differential. 
This strategy has often been adopted, beginning with Einstein fl!J]) and including 
many of the Hamiltonian formulations .IAI21USMI But such a modification is necessar¬ 
ily non-covariant. For our formalism we want to keep our Lagrangian covariant. 
Furthermore a Lagrangian boundary term would modify the boundary condition 
on the whole 3-dimensional boundary of the spacetime region, thus inducing the 
same type of modification on the spatial boundary at large spatial distances as on 
the initial time hypersurface. However we want the freedom to adjust the boundary 
condition on the 2-dimensional boundary of the spacelike region <9E independently 
of the type of initial conditions imposed within the initial time hvpersurface E t . 
Thus for our objectives we turn to the Hamiltonian boundary term|f| 

Note that the Hamiltonian m has two distinct parts; each plays a distinct 
role. The proper density Z^H^, although it has vanishing value on shell, generates 
the equations of motion, whereas the boundary term B{Z) determines not only the 


s In the end it turns out that our favored Hamiltonian boundary term for GR is related to one 
induced by a Lagrangian boundary termE^ 
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quasi-local value (1113D but also the boundary condition. Now we should make note 
of a very important fact: the boundary term can be adjusted— without changing 
the Hamilton equations or the conservation property (11091) . Thus one can replace 
the 2-form B(Z) = izp A p inherited from the Lagrangian by another. 

Such an adjustment is in one respect just a special case of the conserved Noether 
current ambiguity (i.e., for any 2-form x, J and J' := J + d\ are both conserved 
currents (3-forms) if dJ = 0, even though they define different conserved values). 

However here, in this Hamiltonian case, any such adjustment modifies— in par¬ 
allel —not only the value of the quasi-local quantities but also the spatial boundary 
conditions. Thus the boundary term ambiguity is under physical control: each dis¬ 
tinct choice of the quasi-local expression given by the Hamiltonian boundary term 
is associated with a physically distinct boundary condition. 

In order to accommodate suitable boundary conditions we found that, in general, 
one needs to introduce on the boundary for each of the dynamical fields certain 
reference values p, (p, which represent the ground state of the held—the “vacuum” 
(or background Held) values. This is necessary in particular for Helds whose natural 
ground state is non-vanishing; the spacetime metric is such a field. 

We take our boundary terms to be linear in A tp := p — (p, A p := p — p, so that 
they (and thus all the quasi-local quantities) vanish if the fields take on the ground 
state (reference) valuesQ We presume that the reference values (like Z ) are not 
varied: dp = 0 and dp = 0, consequently dAp = dp, dAp = dp. 


14.4. Covariant-symplectic Hamiltonian boundary terms 

To find an improved Hamiltonian boundary term for (II 101) first drop the one inher¬ 
ited from the Lagrangian. Examine the boundary term generated in the variation 
of the 3-form part of the Hamiltonian (II 101) : it is —izpAdp + qdpAizP- This invites 
us to add a suitable complimentary boundary term. In this way we were led to the 
boundary termJ^SEHH! 


B(Z) :=i z { t 


A A p — qAp A iz 



( 120 ) 


Then the associated variational Hamiltonian boundary term becomes 


dU{Z) ~ d 


izdp A A p 1 J —A p A izdp 

—izApAdp) \ dipAizAp 


( 121 ) 


Here for each bracket independently one may choose either the upper or lower term, 
which represent essentially a choice of Dirichlet (fixed field) or Neumann (fixed mo¬ 
mentum) boundary conditions for the space and time parts of the fields separately!)] 


t Some authors use the terminology regularize. 

u There are more complicated po ssib ilities, “mixed” choices involving some linear combination of 
the upper and lower expressionsWe do not have any specific physical examples, but mixed 
boundary conditions may be of interest in certain cases. 
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Iii each of these cases the boundary term in the Hamiltonian variation has a 
certain symplectic structure which pairs certain control quantities— i.e., the inde¬ 
pendent variables—-with certain associated response quantities—the dependent vari¬ 
ables. (For discussions of this paradigm of the symplectic structure associated with 
variational principles which we have found to be illuminating see Refs. iH). 
The symmetry of the above expressions under an interchange of “control” and “re¬ 
sponse”, formally S —> A, A —> —6, is noteworthy. 

Thus, although it is not so well known, when the issue is examined one can 
readily see that there are many choices of boundary conditions and consequently 
really many different expressions for energy in classical field theory. This is true 
especially for gravitating fields. Actually this sort of thing is not unusual in physics; 
in particular one can compare the situation with that in thermodynamics (which has 
several physically meaningful energies: internal, enthalpy, Gibbs and Helmholtz). 

Nevertheless it should be noted that one of our boundary term expressions stands 
out: for any field which allows trivial reference values, <p = 0 = p, one boundary 
term choice vanishes (the lower choice in each bracket). Such fields, with this choice 
of boundary condition, make no explicit contribution to the quasi-local boundary 
term. This particular boundary term has another virtue: for any field with gauge 
freedom it is the only gauge invariant choice. Thus there is a certain preferred 
boundary expression—and thus a preferred boundary condition —for this large class 
of fields, a class which includes all the physical fields of the standard model. There is, 
however, a quite important exception: gravity, more specifically, any gravity theory 
formulated in terms of dynamic spacetime geometry which includes the spacetime 
metric as a dynamical field. The natural reference choice for the metric is not a van¬ 
ishing metric tensor but rather the non-vanishing Minkowski metric. Consequently 
one must have, in general, a non-vanishing Hamiltonian boundary term. 

15. Standard asymptotics 

This section is concerned with suitable asymptotic conditions for our classical fields 
at spatial and null infinity. It also includes a discussion of energy flux. 

For spatial infinity, the issue of asymptotic conditions was first investigated in 
GR by Regge and Teitelboimp 37 with later refinements by Beig and O MurchadhaP 
and then Szabados i 132 ! 13 ^ We have developed a similar idea for general fields. 

15.1. Spatial infinity 

For finite regions these boundary terms in the variation of the Hamiltonian tell 
us exactly what needs to be held fixed (i.e., “controlled”). For asymptotically flat 
regions, however, one should take into account the asymptotic fall off rates. The 
various boundary terms we have constructed enable the Hamiltonian to be well 
defined on the phase space of fields with suitable asymptotic behavior for all typical 
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physical fields. 

For the fields it is sufficient^ to take the respective asymptotic fall offs for even 
and odd parity terms to be 

Aip ss <D + {l/r) + 0~~(l/r 2 ), Ap « 0”(l/r 2 ) + 0 + (l/r 3 ). (122) 

Parity here means the parity of the components in an asymptotically Cartesian 
reference frame. The 2-surface area element has odd parity, so even parity 2-forms 
automatically have vanishing 2-surface integral. 

For asymptotically fiat spaces the displacement should asymptotically be a 
Minkowski Killing vector, i.e., an infinitesimal Poincare displacement. It is suf¬ 
ficient to take 


Z* w Z% + XqvX 1 ' + 0+{l/r) + 0~ (1), (123) 

where (in terms of asymptotically Minkowski coordinates) Zq is a constant trans¬ 
lation parameter and Xq V = Aq 11 ^ is a constant asymptotic infinitesimal Lorentz 
boost/rotation parameter. 

With the asymptotics (112211123D it is straightforward to check that for any of 
the boundary term choices (I120D all of the quasi-local quantities have finite val¬ 
ues; furthermore, for any of the choices all of our Hamiltonians are differentiable 
on the specified phase space—since the respective boundary terms in the varia¬ 
tions of the Hamiltonians (11211) vanish asymptotically. Thus our Hamiltonians are 
generally well-defined on a large phase space which includes physically interesting 
solutions. At asymptotically-flat spatial infinity the aforementioned asymptotics 
are physically reasonable. Our considerations naturally straightforwardly extend to 
asymptotically anti-de Sitter spaces. Here the details are omitted. 


15.2. Null infinity 

Let us next consider what can be expected if the boundary of our 2-surface dT, 
approaches future null infinity. Long range radiation fields (e.g., electromagnetism) 
have slower fall offs, like Ap « dp « 0(l/r). Then the boundary terms in the 
variation of our various Hamiltonians will not vanish, so the Hamiltonian is no 
longer functionally differentiable. This seeming calamity is actually providential -it 
is directing us to additional physics contained within the formalism, namely energy 
flux expressions (HHHHI] 


v Sufficient, but not necessary. When examined in detail it can be seen that one really only needs 
conditions on certain combinations of the components, but it is not in the spirit of our treatment 
to break fields up into, e.g., components parallel to and perpendicular to some specific boundary 
surface, etc. Here we are satisfied with a formalism that includes a large class of fields. We leave to 
more specific investigations finding the largest acceptable phase space with the weakest conditions. 
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15.3. Energy flux 

For the flux of “energy” there is a special way of calculating- the analog of the 
classical mechanics calculation (for conservative Hamiltonian systems) of 

5H = q k Sp k -p k Sq k => E:=H = 0 (124) 


under the replacement 6 d/dt, where the remarkable cancelation is a consequence 
of the particular (symplectic) form of the Hamiltonian variation. In the present case 
from (11161) under the replacement <5 —► £z the same type of symplectic calculation 
occurs, and we are left with the respective contributions from our various boundary 
term choices (11211) 


£ z H{Z) = d 


iz£z£ A A p I f -A ip A iz£zP 

—i Z A<P A £zP / 1 £zp£iz&P 


(125) 


(We are presuming that £zf> , £zP vanish.) 

In particular, as we mentioned earlier, for all fields with vanishing reference, 
there is a standard choice of Hamiltonian boundary term, namely the one that 
vanishes. The corresponding energy flux expression is 


£ z 'H{Z) = d {-izp A £zP + s£zp A i z p) ■ (126) 


16. Application to Electromagnetism 

To illustrate these ideas in a familiar setting, we briefly consider vacuum electro¬ 
magnetism in Minkowski space (for the complete details see Ref. 

For electromagnetism in Minkowski space the formalism developed above, with 
the important exception of the “on shell” vanishing of TL^, is still applicable. A 
first order Lagrangian 4-form for the (source free) Maxwell one-form (the U(l) 
gauge potential) A is 

£^=dA/\V-^Z 0 *P/\V, (127) 

which yields the pair of first order equations 

dV = 0, dA - Z 0 * V = 0. (128) 

These are just the vacuum Maxwell equations with *V = Zf 1 F := Z(f 1 dA- 1 hence 
T = —Zq 1 * F, and d* F = 0. (Here Zq is the vacuum impedance, which has the 
value go = £q 1 in our relativistic units in which c = 1. With our conventions, our 
conjugate momentum held V turns out to be the negative of H which was introduced 
earlier in (15fjl) .') With Z = <9 t and the decomposition A = (—</>, A k ) we find that 
izF = izdA = £zA — dizA corresponds to F ok = A k + d k (f = —E k . The magnetic 
held strength is F\j := diAj — djAi =: eij k B k . Hence Voi = —Zf 1 * F oi = —p(f 1 B i , 
and Vij = — Aq * F i:j = -~£ 0 eij k E k . The natural reference choice is A = 0 = V. 
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The Hamiltonian 3-form is 

U EM (Z) = -i z AdV-dAM z V + i z (^Z 0 *V A V^J + dB EM . (129) 

In the usual tensor index notation the volume density part has the form 

n EM = $d k n k + ±(diAj - djAi)e iik H k + ^-ttS r k - ^H k H k , (130) 

where the momentum conjugate to the 3-vector potential is given the name n k (it 
works out to have in the usual terminology the value —£q E k , i.e., —D k ). By varying 
H k one obtains H k = ^e kl i (diAj — djAi ) = B k , a 2nd class constraint that could 
be used to eliminate the magnetic field, then the Hamiltonian volume density would 
correspond to the familiar energy density ^(eqE 2 + fi Q B 2 ) plus a gauge generating 
term, </><9fc7r fc , which vanishes “on shell”; the scalar potential in this term acts as a 
Lagrange multiplier to enforce the (first class) Gauss constraint d k D k = 0. 

Let us just consider two boundary term choices, namely our preferred choice with 
vanishing boundary term, and the above Hamiltonian 3-form with the boundary 
term 

B EM = i z AV= -<j>n k dS k . (131) 

These two are actually both well known physically, the former corresponds to the 
energy density from the gauge invariant energy-momentum tensor (|8|) . and the 
latter is the energy density of the electromagnetic canonical energy-momentum 
tensor ©• Here our interest is not in the field equations but in the total differential 
term which, upon integration, becomes a boundary term indicating the boundary 
condition. Briefly: for the choice with vanishing Hamiltonian boundary term, the 
total derivative term in the variation of the Hamiltonian is 

- d(i z ASV + 6A A i z V ) ~ d k {(/)5ir k - e kij 6AiHj), (132) 

which tells us that one should hold fixed on the boundary of the dynamical region 
the normal component of the electric field and the surface parallel components of 
the vector potential (the gauge independent part of which determines the normal 
component of the magnetic field). On the other hand, for the Hamiltonian including 
the boundary term (11311) , one finds that the total differential in the variation of the 
Hamiltonian is now 

d(i z 5AV -5AA i z V) ~ -d k {5^ k + e^SAiHj). (133) 

This is the same boundary condition in the vector potential/magnetic sector but 
now for the electric sector one should instead hold fixed on the boundary of the 
region the scalar potential. 

The physical meaning of such boundary conditions are well known. Fixing the 
normal component of the electric field on the boundary corresponds to fixing the 
surface charge density. An instructive physical example is a parallel plate capacitor. 
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One can use a battery to charge up a capacitor with a moveable dielectric. Dis¬ 
connect the battery and measure the work needed to remove/insert the dielectric 
(the potential varies but the charge is fixed, no current or power flow). Alterna¬ 
tively leave the battery connected and measure the work needed to displace the 
dielectric—now the potential is fixed but the charge varies, so current and hence 
power flows. The respective boundary terms in the variation of the Hamiltonian are 
(j>8TT k dSk and —S<f>Tr k dSk- Both boundary condition choices are physically meaning¬ 
ful corresponding to real situations. 

Nevertheless for electrodynamics one expression stands out, the one with van¬ 
ishing boundary term. This choice is the only one in which the value of the Hamil¬ 
tonian is gauge invariant. Moreover, this is the only non-negative Hamiltonian 
density. Consequently the associated energy has a lower bound and the system has 
a natural vacuum or ground state: zero energy for vanishing fields. The value of 
the Hamiltonian with this boundary term can be interpreted as the internal energy , 
whereas the other expressions can be regarded as including some additional energy 
on the boundary of the system associated with maintaining the boundary condi¬ 
tion. The associated electromagnetic energy flux expression from our formula (11261) 
reduces to just the usual Poynting energy flux: 

£zU em = d \—izA A (diz + izd)'P — (diz + izd)A A izV] 

= -d(i z F \i z V) =d(-E i H J dx i A dx j ) . (134) 

Clearly this choice, associated with fixing the normal components of the electric and 
magnetic fields on the boundary, is preferred; it is the one suitable for most physical 
applications. It gives the usual energy density and Poynting energy flux. Similarly, 
for all other fields— except for dynamic spacetime geometry —there is available a 
standard Hamiltonian (the one with vanishing boundary term contribution ) associ¬ 
ated with a certain preferred boundary condition. 

17. Geometry: covariant differential formulation 

In the discussion of our covariant Hamiltonian approach, up to this point (except 
as was specified for a couple of specific examples), there has been no need to make 
any restriction on the type of geometry for our manifold. Here in this section we 
discuss the specific sort of dynamic spacetime geometry that we will consider and 
relate it to the gauge theory paradigm. 

The covariant Hamiltonian formulation can apply to general theories of dy¬ 
namical geometry. Standard references for differential geometry are Kobayashi & 
NomizrP 3 and SpivakP^H 

17.1. Metric and connection 

For the dynamical spacetimes that we consider there are two basic geometric ideas: 
a metric tensor g = g^d^ ® d v (which determines length and angle), and paral- 
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lei , associated with parallel transport , covariant derivative and connection , i.e., we 
consider metric-affine geometry. 

The metric gives the causal structure, arc length, area and volume. Furthermore 
it is used to raise and lower indexes (i.e., it determines a specific isomorphism 
between tangent and cotangent vectors). It also provides the paths of extremal 
length (geodesics). For the 4D spacetimes of interest here, the metric has the 
Lorentz signature. Given such a metric there is naturally defined an associated 
symmetry group, the group of local Lorentz transformations: L £ SO(l,3) => 
g(LX,LY)=g(X,Y). 

For the other structure, let e M for /i = 0,1, 2, 3 be a basis for spacetime vector 
fields. The covariant differential V of each basis vector is a vector valued one-form, 
hence some linear combination of the ep’s with one-form coefficients: 

vep = e a r a p, r a p = r a p i dx i , ( 135 ) 

called the connection one-forms. 

The covariant differential of a vector field V = V M is 

W = V(e M V*) = (Ve M )^ + = e v T\V^ + =: (136) 

Its components are determined by the operator D : 

DV> i :=dV tl + r^Ar, (137) 

which extends, as indicated, to vector valued forms. 

The notation automatically antisymmetrizes: 

V 2 V = V 2 ( e/i l/ M ) = V( e/i W M ) = e^D 2 V li 

= e M [didV^ + A V v ) + r^ A A (dV x + T x a A V a )] 

= e M (dr^ + r^ A A r A „) AF = e^ v A V v t (138) 

where 

R^ v := dT"* + r^A a r A „ = \^R IJ vl] d:C A dx'i. (139) 

is the curvature 2-form 

Exterior covariant differential form notation treats some, but not all, indices as 
differential forms. Rather than work with the operator V on geometric objects we 
often find it more convenient to work with D on their coefficients. The covariant 
differential D can be extended to operate on a tensor valued form of any type, e.g., 

DP a p = dP a p + Y a 1 AP^p~Y^pAP a 1 , (140) 

D 2 P a l6 = P“ 7 A P 1 p - PC p A P“ 7 . (141) 

Generically, arranging the components of a tensor valued form p as a row vector 
we get 


P 

a , 


Dip = dip + T “p A pa, 
D 2 p = R a p A pCTa ^, 


(142) 

(143) 
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for some appropriate representation matrix o a ^. 

The special case of the vector whose components are the coframe one-forms 
fla _ e <x ,fl x i yjgijjg the torsion 2 -form: 

T a := Dd a := dd a + Y a p A = ^T a ij dx i A dx j . (144) 

On the coframe i?“, D 2 gives an important special case of (11431) : 

DT a = D 2 d a = R a f 3 A ftP, (145) 

which is known as the first Bianchi identity. 

Applying D to the curvature 2-form R a p gives 

DR a p := dR a p + T “ 7 A R?p - Y 1 p A R a 1 

= d(dY a p + T Q 7 A r/ p ) + T “ 7 A (dTp + Y^ a A Y a p) 

-Y^p A (dT “ 7 + Y a a A r CT 7 ) = 0, (146) 

by explicit calculation. This is the second Bianchi identity. 

With the metric tensor, Dg^ defines the non-metricity 1-form. 

Dg, lv \= dg^ v - Y x /J ,g Xu - Y x v g llX - (147) 

Correspondingly, we have 

D 2 g^u = -R x ^g\ v - R x ^gfj.\- (148) 

17.2. Riemann-Cartan geometry 

Here we are interested in particular in the special case where the geometry can 
be regarded as a local gauge theory of an appropriate spacetime symmetry group. 
With due consideration given to the understanding of both the geometry of and 
physics in Minkowski spacetime, the appropriate choice for the symmetry group 
is the inhomogeneous Lorentz group, generally referred to as the Poincare group. 
This is both the symmetry group for the spacetime of special relativity and the 
group used to classify elementary particles in terms of mass and spin. The group 
is a semidirect product of the translation group and the group of rotations/Lorentz 
boosts. The Noether conserved quantities associated with these global symmetries 
are energy-momentum and angular momentum/center-of-mass momentum. The 
type of spacetime geometry with local Poincare symmetry is known as Riemann- 
Cartan geometry. 

In Riemann-Cartan geometry the connection is assumed to be a priori metric 
compatible , Dg^ = 0, via (11481) this gives R a p = R[ a p] (he., a Lorentz Lie alge¬ 
bra valued two-form). For our purposes it is convenient to use the orthonormal 
frame gauge condition, then the metric components are constant and dg^ = 0 , 
so, via (11481) . this gives Y a p = P[ Qi s] (be., a Lorentz Lie algebra valued one-form). 
The geometry has in general nonvanishing torsion and curvature. The metric infor¬ 
mation is encoded in the orthonormal coframe which has local Lorentz gauge 
freedom. 
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Riemannian geometry is a special case with vanishing torsion, T a = 0; such 
a connection is called symmetric. Then the geometry is given by the curvature, 
which is generally non-vanishing, so the parallel transport is path dependent. On 
the other hand another special case is teleparallel geometry, which has a vanishing 
curvature 2-form, R a p = 0. This is referred to as flat. Parallel transport is then 
path independent, nevertheless it is generally nontrivial -being dependent on the 
torsion, which is generally non-vanishing. 

17.3. Regarding geometry and gauge 

Note the respective similarities in the form of the commutator of the gauge covariant 
derivatives for the flat space U(l) phase case (V M </> = d^cj) + ieA^cj)) and the Yang- 
Mills case (V M V> = d^ip + iqA^Tptjj) compared with the spacetime geometric case: 

[V^, V„] <j> = ieF^cj), (149) 

(150) 

[V„, V,] V a = R a p^VP - TVV/“. (151) 

Here is the covariant derivative ; the latter relation is called the Ricci identity. 
On the right hand side the respective gauge field strengths appear. One can see 
that for spacetime the curvature is the Lorentz field strength, and the torsion is the 
spacetime “translational” field strength, associated with the generator of infinitesi¬ 
mal translations, the directional derivative. 

Riemann-Cartan geometry is ideally suited to admit an interpretation as a local 
gauge theory of the symmetry group of Minkowski space, the Poincare group. (In 
the standard Riemannian GR formulations the torsion is a priori assumed to vanish, 
then gravity does not look much like a local spacetime symmetry gauge theory. 
Teleparallel geometry can be regarded as a gauge theory for translations.) 

We see that, when suitably formulated, gravity has both Lorentz/rotational and 
translational “vector potentials” which are similar to those of the Maxwell/Yang- 
Mills theories. 

17.4. On the affine connection and gauge theory 

The “connection” one-forms for “translations” and “Lorentz” transformations can 
be packaged together in a way that offers some further insight into their essential 
similarities and differences. 

The Poincare transformations on Minkowski spacetime 

V a '=A a 'pV p + A a ' (152) 

can be conveniently represented in matrix form as 


(153) 










July 28, 2015 0:26 


World Scientific Review Volume - 9.75in x 6.5in 


100GR'Ch4'150726vl page 44 


44 


C.-M. Chen, J. M. Nester, R.-S. Tung 


Then the matrix product 

( Ai A\ \ / A 2 A 2 \ _ / AiA 2 AiA 2 + Ai \ 

Vo i J\ o i J ~V o i ) 


(154) 


reflects the semi-direct product structure. This matrix representation for infinites¬ 
imal Poincare transformations has the Lie algebra 

h ai \ / h A 
o o J ’ Vo 0 J 


[h,h] ha,2 — hai 
0 0 


(155) 


Now a connection can be viewed as a Lie algebra valued one-form. The spacetime 
“translation” and “Lorentz” connections can thus be neatly packaged in terms of 
the above Poincare Lie algebra matrix representation: 


u j := 


r i? 
o o 


(156) 


The associated “curvature” Lie algebra valued 2-form 


fl:=dw + «Aw 


/ dr + r a r 

V o 


dd + TAd\ 

0 J 




(157) 


includes the spacetime curvature and torsion 2-forms in one package. Furthermore, 
the Bianchi identity for this Poincare Lie algebra curvature matrix, 


0 = M1 :=(ffl+wAfl-OAw (158) 

dR + T/\R-RAT dT+ T AT - R Rtf\ _ ( DR DT-RAd 

0 0 ) ~ V 0 0 

unifies the first and second spacetime Bianchi identities. 

This packaging shows similarities between the connection and coframe one-forms 
and the curvature and torsion 2-forms, but also some clear differences inherited 
from the semi-direct product structure of the Poincare group. The gauge theories 
of Yang-Millfi^ and Utiyam d 139 * 140 ! also have the fl = dui + ui A u and DCl = 0 
form, but the groups do not have a semi-direct product structure. 

Although we find this formulation quite helpful for seeing how the coframe plays 
the role of the “vector potential for translations”, we will not use it below in our 
treatment of the PG. For our purposes we consider local Poincare transformations 
to be Lorentz transformations of the coframe plus local spacetime diffeomorphisms. 



18. Variational principles for dynamic spacetime geometry 

In this section we develop the second order variational principle for gravitating 
material and internal gauge fields along with their associated Noether currents and 
differential identities. The spacetime is assumed to have Riemann-Cartan geometry, 
i.e., we are considering the Poincare gauge theory of gravity (PG). 
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We are considering geometric gravity that can be regarded as a gauge theory 
for the Poincare group. Several authors have considered such theories, see, e.g., 
Refs. I7l8l48l53l57l811 . 

We wish to consider the conserved Noether currents and differential identities as 
well as the field equations for dynamic spacetime geometry and gauge interactions. 
Here in this section we first work with the usual 2nd order type Lagrangian, since 
for that case certain expressions take a simpler form and the arguments are more 
transparent. Having established these results we can then present more briefly 
the analogous first order version which is the basis for our covariant Hamiltonian 
expressions. 


18.1. The Lagrangian and its variation 

Rather than begin with the Lagrangian 4-form of the typcE 


C = £(p, iT, T“s, A p ; dr“/ 3 , dA p ), 


(159) 


we take the more covariant form 


C = £(p, T a p ,A p -Dp, T p , R a p, F p ), 


(160) 


which is no less general. Here p is a generic /-form source field with total covariant 
differential (the factor ordering is suitable for a matrix representation with p as a 
row “vector”) 


Dp = dp + T a p A po a P + A p A pT p . (161) 

The torsion 2-form, curvature 2-form and gauge field strength were given ear¬ 
lier (1144111391 EH). The respective variations are 


SDp = DSp + 8T a p A paj + SA P A p T p , (162) 

ST p = DSd p + <SF*„ A rT, (163) 

5R a p = d6T a p + ST a x A T x p + T a x A 4T a ^ = D5T a p, (164) 

SF p = dSA p + C p qr A q A 5A r = D6A p . (165) 


w This is sufficiently general to include all the fundamental fields of the standard model, with 
including the Higgs and the fermions and A p the U{ 1) x SU( 2) x SU( 3) gauge vector potential. 
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The variation of the total Lagrangian 4-form (11601) is: 


SC = SDp A 


dC 

dDip 


dC dC 

ST P A + SR a p A + SF P A 


dTh 


+< *’ A | + w ' ,A ^ + 4r “-> A ar 


dR° 

dC 


dC 

8Fp 


SA p A 


dC 

8Ap 


BC dC 

= (DSip + ST a a A poj 3 + SA P A ipT v ) A —— -b Sip A — 

‘‘dip 


+ (DSd p + ST P V A r) A 


<9£ 
d Tp 


5d p A 


acv 

oc_ 

dd p 


+ DSTa P^^-+ Sra ^ ^ 


+DSA P A 


5J?° 

<9£ 


ar° 


A 


dC 


dFP ' '' <9Ap' 

“Integrating by parts” and rearranging gives 
dC 


SC = D \Sip A 


siv 


dC 

Sd p A ^7 + tfT o 0 A 


A ( —<;D 


dC 


dTP 
dC 


dDtp dp 


+ Sd p A ID 


dC 

dR a c 

dC 


+ 5 T a p A ( D 
+SA p A ( D 


dC 


dC 


dR< 

dC 


dv 
dC 


+ (pcr a p A 


dTP 

dC 


- SA p A 

dFP J 

dC \ 
di!)P ) 


+ T~: -b (pTp A 


dDip 
dC \ 


■ A 


<9£ 

5T“ 


(166) 


(167) 


(168) 


' dAP r p '' dDip J ’ 

which yields the conjugate momenta and Euler-Lagrange variational derivatives 
according to the pattern: 

SC = d(Sp A p + Sd a Ar„ + ST a p A pj 3 + SA P A V p ) 


SC SC SC 

+ 5 ip/\- — b A—— \-SF p A—— 

r - s§° p sr a f 


SAP A 


SC 


(169) 


Stp Sd a p ST a p SAp 

This is the key variational relation. According to Hamilton’s principle, the 2nd 
order field equations are the vanishing of the Euler-Lagrange expressions named 
in (11691) and explicitly displayed in (11681) . 


18.2. Local gauge symmetries, Noether currents and differential 
identities 

The Lagrangian 4-form C (11601) should be “gauge” invariant under the local space- 


time gauge transformations: 

Aip = l a p ipoa 13 + a p ipT p - £ Z ip, (170) 

AiT = - £ z d p , (171) 

AT a p = -Dl a p- £ z T a p, (172) 

A A p = —Da p - £ z A p , (173) 
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where the 6 l a p control an infinitesimal Lorentz rotation of the space-time frame 
d a , the a p control an “internal” gauge and Z is a spacetime vector field which 
determines a “local translation”. The Lie derivative £z is given by di z + i Z d on the 
components of our fields. This action on the components and on the basis one-forms 
correctly represents the Lie derivative on geometric objects. Under (11701) (11731) the 
Lagrangian 4-form C (I160D . if it depends on position only through the indicated 
fields, should change according to 

A C = -£ Z C = -di z C, (174) 


which happens to be a total differential because £ is a 4-form. With the special 
variations given by flnnp cup eq. (11691) should be satisfied identically. 

One may collect all the total differential terms on the l.h.s. of the identity giving: 

dJil-f, c>, Z) s A A + Ail» A A + if", A A- + AAP A A (175) 


where 

J(Fp,a p ,Z) 


—iz£ — Ay> A 


dC 

dDip 


- AiFA 


dC 

dT p 


AF^A 


dC 

dR a p 


- AA p A 


dC 

dF pl 


(176) 

is the generalized total current 3-form. Using (11701) (11731) we have in more detail 
(recall the canonical stress tensor (1551) and E := Q£q — L) 

BC 

J{l a p , Z) := -i z C + {£zv ~ l a p^paj - a p ^T p ) A — 

3C r)C 

+(£z# a - rptf) a ^ + {£ Z T a p + Dl a p) A — 
d£ 

+(£ z A p + Da p ) A—. (177) 


The second Noether theorem differential identities may be obtained from (11751) 
by comparing the coefficients of a p ,l a p, Z p ,da p ,dl a p,dZ p on both sides. The re¬ 
sults will be covariant, but the computation will not be manifestly so—because of 
the Lie derivative terms. However, it so happens that the Lie derivative differs 
from a covariant operation only by a gauge transformation, as the following short 
computations reveal: 


£ z y := di z p + izdip = Di z <p + izDip - T a p(Z)(pa a p - A p (Z)ipT p: (178) 
£ z ^ := di z 0" + i z dd p = Di z d p + i z D- T p v {Z)r, (179) 

£ z T a p := di z T a p + izdT a p = i z R a p + D(T a p(Z )), (180) 

£ z A p := di z A p + i z dA p = i z F p + D(A P (Z)), (181) 


(in the last two relations on the r.h.s., the D is formally defined by treating A P (Z) 
and T a p(Z) as tensors). Thus the translation vector field Z induces a modification 
to the gauge parameters 

l' a p:=l a p+T a p(Z), a ,p :=a p + A p (Z), 


(182) 
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which effectively replaces the non-covariant £ z by the “covariant Lie derivative” 
defined by 


Lz := Diz + izD, 


(183) 


on the “normal” fields and by 


L z T a p := izR a p, L Z AP := i z F p , (184) 


on the connection one-forms. 

The gauge transformations (I170D (1 1731) then take the manifestly covariant form: 

A ip = l ,a pipa a P +a ,p (pT p -L z <p = l la pipaa 13 +a ,p ipT p - Diz^-izDip, (185) 
A0" = l'^ v r + 0 - Lz#* = + 0 - DZ p - i z D-d p , (186) 

A T a p = - Dl' a p + 0 - L z T a p = -Dl' a p + 0 - i z R a p , (187) 

A A p = 0 - Da' p - L z A p = 0 - Da' p - i z F p . (188) 

If one specializes to matter fields which are not form fields—which is the only kind 
of matter that we know of physically—then one can see here a striking pattern which 
supports our identification of the geometric gauge fields. From the r.h.s. expressions 
one finds that most of the fields transform algebraically under gauge infinitesimal 
internal gauge transformations a p ; the only field that transforms with Da p is the 
internal connection one-form (a.k.a. gauge vector potential) A p . Similarly, most of 
the fields transform algebraically in l' a p; the only field which transforms according 
to the differential Dl' a p is the spacetime connection one-form V a p. Moreover (in 
the “physical” 0-form matter case) the only field having a differential- DZ P — 
rather than an algebraic transformation formula under the infinitesimal spacetime 
displacement Z p is the coframe one-form This is one more way of seeing that 
the coframe one-form can be identified as the “translational gauge vector potential”. 

With the above reparameterization the generalized total current 3-form takes 
the form: 


J(l' a p, a' p , Z) := -i z C + (Lz<P ~ l' a p<pi rJ - a*<pT p ) A — 

r)C r)C 

+(Lz <?“ - l ,a p&) A — + {L z T a p + Dl' a p) A — 
dC 

+{L z A p + Da' p ) A—. (189) 


To extract the differential identities from (11751) it is best to write J ( l' a p,a ' p , Z) 
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as terms algebraically linear in l' a p,a' p ,Z plus a total differential: 






(190) 


The total differential will later be related to the total energy-momentum. For now 
merely note that it does not contribute to the l.h.s. of (117511 as d 2 = 0. 

Internal gauge symmetry. With Z = 0 = l' a p the general result (11901) reduces to 
the expression we considered earlier for internal gauge symmetry of the Yang-Mills 
type HMD, an d we again obtain just as before (recall the argument leading to (1761) 
and (1701) 1: 



(191) 


Local Lorentz gauge symmetry. In the same fashion, with Z = 0 = a' p we obtain 


from (11751) and (11851) (11871) 



Equating the coefficients of l ,a p and Dl' a p (keeping in mind that l' a p is antisym¬ 
metric) leads to the algebraic and differential identities: 



(193) 


(194) 


Formally these two relations are quite similar to those found for the internal sym¬ 
metries (11911) . Consequently local Lorentz symmetry is in certain respects rather 
like an internal gauge symmetry. 

The conditions dC/dA p = 0 = dC/dT a p mean (as we expected) that how the 
internal and Lorentz gauge potentials can appear is quite restricted—i.e., only via 
the covariant derivative or the associated field strength. With these conditions we 
note that the generalized current has the neat form 



(195) 
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Local diffeomorphism invariance. This was already considered for general theo¬ 
ries in first order form in section rid.il Here we consider local translations in second 
order form while distinguishing between the source, gauge and geometric variables. 
This case again has some similarities to the internal and Lorentz symmetries but 
also some striking differences. With l' a p = 0 = a' p in (11751) we have 


dj(0, 0, Z p ) 


d(Z p J^ = DZ P A Ju + Z P DJ^ 

r r r r 

-(Di z <p + izD<p) A -(. Di z ^ + izDd p ) A —- 


-izR a t 3 A 


SC 


izF p A 


SC 
JaP' 


(196) 


The coefficient of Z p on both sides of (I196D is the differential identity associated 
with translation invariance (related to energy-momentum) while the coefficient of 
DZ P provides a new algebraic expression for J jL in terms of the Euler-Lagrange 
variations: 


_ _ SC SC 

^ A 6<p Sd p ■ 


(197) 


For ordinary (i.e., 0-form valued) matter fields the first term vanishes, leading to 
the especially simple and intuitively reasonable relation: 




6C 

6dh ' 


(198) 


We emphasize that such a neat formula for the translational current 3-form is only 
possible is one regards the coframe as a dynamical variable. It is also noteworthy 
that this remarkable relation—which does not restrict how the translation gauge 
potential appears in the action—was obtained from the coefficient of DZ P , whereas 
the corresponding relation in both the internal and local Lorentz cases put quite 
severe limits, viz. (iTsm mm on how those gauge potentials could appear. 

The identity (11971) permits J(l' a p , a ' p , Z) (11951) to be written in the nice sym¬ 
metrical form 


J(l' a p, a 


tp ry\ ■ . SC SC /a SC 

p Z ) s _ wA _- !z ^__, 

+d (iz<fi A ttft- + : + l' a 


v'P 


JT- ( 199 ) 

5A p x ’ 


dDip 


DTP 


dC 

dR a f 


v'P 


dC \ 
dFP J ' 


Restricting to the case where the source fields are 0-forms, we have 


J(l' a p,a ,p ,Z) = -Z p ^L-l ,a 


SC 

SdP 


SC 

ST^f. 


— a 


/p 


SC 

~SAp 


( 200 ) 


+d (z^ + r- dc 


CP 


\ dTP 


d R° 


dC 

dFP ) 


V 


This result displays the purely gauge nature of the current—it is especially note¬ 
worthy that there is no explicit appearance of the source field or its field equation. 
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Note that if we take the variational derivatives as field equations, the numerical 
value of J will be entirely from the total differential term. When the 3-form J 
is integrated over a 3-dimensional region this total differential becomes, via the 
fundamental boundary theorem, (1341) . an integral over the 2-dimensional boundary 
of the region. In other words, the value is quasi-local. 

Our results here are an application of Noether’s ideas. As expected from her 2nd 
theorem, with a local symmetry the conserved current becomes a differential iden¬ 
tity. Furthermore, we also displayed here detailed results that exactly reflect her 
remarks about verifying and generalizing Hilbert’s assertion regarding the lack of a 
proper energy-momentum. Our generalized current expressions (11991) and (12001) are 
linear combinations of Euler-Lagrange expressions plus a total differential. They 
have the usual conserved current ambiguity: the total differential does not con¬ 
tribute to the conservation law, it can be adjusted without affecting the conserva¬ 
tion property, nevertheless it affects the value of the associated conserved quantity. 
The second order Lagrangian formalism has no way to control this ambiguity. 


18 . 3 . Interpretation of the differential identities 

Let us now specialize and consider the customary “minimal coupled” decomposition: 

Aotai = Co r(., ■T ,.,.;., T", R a p ,.) + C A (., .,.;.,.,., F p ) 

+£ v ( V) 0**,.,.; Dp,(201) 

Each of these separate Lagrangian pieces is a scalar valued 4-form. So the Noether 
identities we have obtained can be applied to each piece. However, it must be kept 
in mind that for the separate pieces our variational derivatives are, in general, no 
longer field equations. There is, however, one exception: 


6Cm total 


( 202 ) 


6<p 5<p 

which vanishes “on shell”. It should also be noted that for each separate piece of 
^totai many of the terms in the various identities dun m ma sszp will vanish 
trivially. Consider (11911) : for C$r it is trivial; for C A and C v it reduces to the results 
obtained earlier ( 17511 . 

Let us introduce some suitable names for certain expressions. Specifically for 
the material source the energy-momentum and spin density 3-forms are defined 
respectively by 


dC v 


{-Pa o 

V ‘ dT 


ot.fi 


with an analogous expression defining T^. 

Then (11941) for C v and C A becomes, respectively, 

a r 

D&°f + iT A - 00 A T£ = 2po afi A —= 0 (on shell), 

^ ^ ^ Oif 

2 = o. 


(203) 


(204) 

(205) 
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The physical interpretation of these concerns angular momentum conservation (or 
more precisely the exchange of angular momentum with the gravitational field). 
The first term in (12041) is the (covariant) divergence of the source spin density, 
the next two terms (the anti-symmetric part of the energy-momentum density) 
describe the change in “orbital” angular momentum. The second relation is an 
identity which shows that gauge fields have symmetric energy-momentum densities 
and trivial angular momentum conservation relations. For C$y (11941) reduces to the 
identity 




A^-Rfl, 


dv 




dC$T 

dT ai 


+ % A [ D 


r 

d T«1 


=0 


(206) 


which is satisfied iff is a local Lorentz scalar. 

Now let us consider the differential identity part of (11961) . For a 0-form material 
source field it takes the form 

D‘l* - i e J v Air- \ie,R af) A &p a = = 0 (on shell). (207) 


When the source field equation is satisfied this relation describes the exchange of 
material energy-momentum with the gravitational held. The analogous expression 
for the gauge held Lagrangian is a little simpler: 

^-^r A ^ + ie/ p Ai)|^E0 (208) 

the interpretation is similar. 

Finally, for the differential identity of (11961) applied to r, after some straight¬ 
forward calculation, we have the identity 


D 


dC$Y 

ddh 


- 


,r a [d 


dT" 


dC$r \ 

dd” ) 


R af3 A D 


dC$T 

dR a P 


- DT a A i e 


dC$Y 
■ dT a 


= 0, 


(209) 

which is satished if C$y is a scalar valued 4-form constructed out of the coframe, 
torsion and curvature. This is the PG identity that plays a role analogous to that 
of the contracted Bianchi identity in GR. 

To obtain these detailed results we used the dehnition of the various Eulcr- 
Lagrange expressions given in (I168I169D . 


19. First order form and the Hamiltonian 

Here for our general PG with generic matter and gauge sources we briefly present the 
first order form along with the associated Noether currents and differential identities 
and then the associated covariant Hamiltonian including our preferred boundary 
term which yields our quasi-local quantities. 
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19.1. First-order Lagrangian and local gauge symmetries 

For certain purposes we find a first order formulation convenient. The first order 
form of our variational principle for geometry and gauge is 

£ lst = Dtp A p + Dd a A r a + R a P A p a p + F p A V p 

-A (ip, d a , T^,A p ;p, r a , p a p,V p ), (210) 

with T a P, R a P, and p a p being a priori antisymmetric. The variation takes the 
pattern 


<hC lst =: d{ 6 ip Ap + 5d a A r a + c5T a/3 A p a0 + SA P A V p ) 
SC lst 


+ 6 <p A 


Sip 


AA lst AA lst 

Sd a A -P- + <5T Q/3 A 


SC lst . 

H-r— A Sp - 

dp 


SC 1 


St„ 


Sd° 


A 6 t° 


SC 1 


Spa/3 


5T a P 
A Spafi 


SA P A 


5C 1 


SAP 


SC lst 

sv v A6Vp ’ 


( 211 ) 


where 


5C lst 

5T a P 


6 C lst „ dA 

— = - iDp -w 

SC lst dA 

Sd a ~ Ta ~ dd^’ 


dA 


Dp a p ~ + Wapp + d[p A r a] , 

SC lst ^ dA 

= DV v - tttt + pT p p, 


SAP 


dAP 


SC lst 

Sp 

sc lst 

Sr a 

sc lst 

Spotfi 

sc lst 

~5Vf 


= Dip 


= Dd a 

= Ra p 


dA 
Up’ 
dA 
dr a ’ 
dA 

dpoip 


= F p — 


dA 

dK 


( 212 ) 

(213) 

(214) 

(215) 


It is instructive to compare these (and subsequent) relations with the corresponding 
ones in the previous section. Here we have twice as many fields and thus twice 
as many Euler-Lagrange expressions, but, on the other hand, the Euler-Lagrange 
expressions are all linear and much simpler. 

The first order formulation is convenient for imposing certain constraints. In 
particular in order to impose one of the conditions 


R a f 3 = 0 teleparallel connection, (216) 

T a = 0 symmetric connection, (217) 

in the first order formalism we need merely take the potential A to be independent 
of the corresponding conjugate momenta. The momentum then functions as a La¬ 
grange multiplier which imposes the constraint. The related “coordinate” equation 
then loses its dynamical significance and instead becomes a relation for determining 
the multiplier. 
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19.2. Generalized Hamiltonian and differential identities 


To obtain the Noether differential identities in this mode we need, in addition 


to (11851) (11881) . the gauge transformations of the conjugate momenta: 

A p = -l ,aP (r a pp - a ,p T p p - L zPl (218) 

At 0 = -l' a pT a - L z rp , (219) 

A pj = I'^pJ - r aP ^ - L ZP J, (220) 

A V q = -V p C p qr a' r - L z V q , (221) 

where L z = Di z + i z D. These results were deduced from relations like 

A(p A Dip) = Ap A Dp +p A A(Dtp), (222) 


using the fact that p A Dtp is a scalar valued 4-form. 

As before AC is a total differential: —di z C. Taking the total differential terms 
to the l.h.s. in (12111) gives 


dH(l' a p,oF,Z ) 

A <S£ lBt . aa SC lst ^ a8 SC lst . sc lst 
= ApA- ^A — + AT^A^ + AA p A — 

SC lst . AC lst . 5C lst . SC lst 

+—— A A p -\——— A At a H— -A A p a0 -\—— A AT’ V 


Sp 8r a 5p a p 

where the first order “generalized current” is 


SV v 


(223) 


H(l' a 0 ,a' p ,Z) 

:= -i z C lst - Atp A p - Ai9 a A r Q - AT a/3 A p a/3 - AA P A (224) 
= izA - i z (Dtp Ap + Dd a A r a + -R a/3 A p a0 + F p A V p ) 

+(L z p - r 0P oJ - a ,? VT p ) Aj) + (F z tf a - /'V^) A r a 
+{L z T af3 + Dl' afi ) A p a0 + ( L z A p + Da' p ) A V p (225) 

= i z A + A - £>$“ A - F a/3 A - F p A 

+(F>i z ^ - l ,a pipaj - a'PpTp) Ap + (Di z 0 a - J' 0 ^) A r a 
+F>Z ' a/3 A p a/3 + F>a' p A Vp (226) 

= izA + S-Dp A izp - Dd a A i z r a - i ? a/3 A i z p a p - F p A i z V p 
+<;i z tp A Dp - i z d a DT a 

—l' a P(Dp a0 + ip<j a p A p + A T a ]) - a p (DV p + (pT p A p) 

+D(i z ip A p + i z d a r a + l' afi p a p + a' p V p ). (227) 


This expression is not just a Noether current, it is, as we already justified quite 
generally in subsection 114.21 the (generalized) Hamiltonian 3-form , i.e., the canon¬ 
ical generator for internal, local Lorentz and local spacetime displacements (which 
includes the time evolution for any choice of time). Nevertheless, like the 2nd order 
current, one still has the various differential identities. 
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Local internal gauge symmetry. Equating coefficients of a' p and Da' p in (I223D 
gives the algebraic and differential identities: 


dA 

dA p 


= 0 , 


t£ lst m SC lst 
— D—r ~.— = ipT q A 


SC 


1st 


6 A* 


Sip 


Sp 


A Tqp- 


SC 


1st 


sv p 


A V r C r 


(228) 


which can be compared with G2D- The signihcance is the same, but now the 
r.h.s. includes also Euler-Lagrange variations w.r.t. momentum variables. If these 
equations are imposed, the expression has the same form as (HSH). 

Local Lorentz symmetry. Equating coefficients of l' a & and Dl ,otl3 gives the alge¬ 
braic and differential identities: 


dA 


= 0, 


dr a P 

sc lst sc lst a sc lst 

-D v = ip(T\n,pi\ A —--1- V[p A 


SC 


1st 


£pa/3 — C a [aP] 

5C ls t 
St^° 


Sip 
A T/3] - 


SC 0 


SC 1 


Sp^ 


A Paj 


Sp 
SC lst 


A ®[aP]P 


Sp 7 


A P/ 37 , 


(229) 


(230) 


which should be compared with (119311 and (119411 . The significance is the same, but 
now the r.h.s. of the latter includes also several Euler-Lagrange variations w.r.t. 
momentum variables. If these relations are imposed the expression has the same 
form as CMD. 

Local translation symmetry. With the decomposition "H(0,0, Z^) = Z^T-L^ + dB , 


Dn{0,0,Z p ) 


= DZ p A Hu + Z^DHft 


= - L z ip A 


AA lst AA lst AA lst AA lst 

- - L Z T^ - L Z A P C 5C 


Sp 


Sd° 


ST a P 


6 A p 


SC lst r 5C lst r SC lst r SC lst r m 

A L z p ——A L Z T a — — — ALzPu/3 — 7 ——A L Z V P . (231) 

Or t) 


Sp 


Sr n 


Spap 


From the coefficient of Z p we obtain a differential identity involving ; this relation 
includes the conservation of energy-momentum. 

We also get a new algebraic identity giving in terms of variational derivatives: 
(compare (11971) 1: 


Up = A- 


SC 1 


SC 1 


Sip 


SdP 


SCist §c ist 


SC 1 


A T afl ~ 

Spap 


A Pa/3 /a 


SC 1 


sv E 


■A'P, 


pm 


(232) 

When the momentum relations are imposed, this reduces to the corresponding ex¬ 
pression for the second order Noether translational current (11971) . Inserting this 
new result into (122711 gives an expression for the Hamiltonian 3-form in terms of 
variational coefficients (which has the same form as (11991) if the momentum relations 
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are imposed) 

H{l ,a 0 ,a' p ,Z p ) 

= -izP A 


SC 


1st 


Sip 


- izd p 


SC 


1st 




-r a0 


sc 


1st 


ST 


— a 


/ P 


SC 


1st 


sc lst sc lst . sc lst 

+<T— 7 — A izv - 7 -A i Z T a - -A %zp a p ~ 

op oTq, Sp a p 

+D(i z p A p + izd a T a + l ,a0 pap + Ot' P Vp). 


SAP 

SC lst 


5V v 


A izVp 


(233) 


This generalized current expression is again an example of applying Noether’s 
analysis. In accord with the second theorem, for local symmetries the current 
conservation expression becomes a differential identity. And again we have a de¬ 
tailed expression that reflects Noether’s remarks regarding a more general version 
of Hilbert’s assertion concerning the absence of a proper energy-momentum conser¬ 
vation law. Moreover there is again the Noether current ambiguity regarding the 
total differential term. However, as explained in section [TO the first order current 
is also the Hamiltonian, the canonical generator of local transformations including 
spacetime displacements. The total differential (boundary) term in the Hamiltonian 
can and should be adjusted, as we have discussed in general terms earlier in sec¬ 
tion [TT31 Then, when varied, the chosen boundary term in the variation variation 
of the Hamiltonian gives the associated boundary conditions. Thereby the usual 
Noether current ambiguity is fixed by the chosen boundary condition. 


19.3. General geometric Hamiltonian boundary terms 


Specializing our general Hamiltonian boundary term expression (11201) to our present 
variables, with the preferred choice for the material and internal gauge fields leads 
to boundary term expressions which explicitly contain only the geometric variables: 

B(Z) = i z j|“ Jat q +A 0“ Aizjj*} 

+ iz^y]ApJ + AT a 0 Ai z ^ P _y^ , (234) 

where the upper or lower line in each bracket is to be selected. A special case of 
this expression, (upper, lower, upper, upper) with f a = 0, was first proposed in 
199lP^ With the above boundary term, the total differential term in S'H(Z) has 
the symplectic form 


f izSd a A At q 1 f A$“ AizSr a j 
\ —izA-d 01 A Sr a J \ —Sd 01 A iz^T a J 

r izST a pAAp a 0 I f AT a pAi z SpJ \ 
+ \ —i z AT a p A 6p a 0 J j ST a 0 A i z Ap a 0 J 


(235) 


The general geometric Hamiltonians evolve gauge dependent quantities, includ¬ 
ing the frame and connection coefficients. Consequently they naturally include 
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terms that are gauge dependent. These terms are in particular those with factors 
of iz£ a However the value of the Hamiltonian boundary term will then include 
a contribution to the energy etc. due to the choice of frame gauge. From such a 
boundary term one could still get the correct physical value, but only if one takes on 
the boundary the particular frame gauge £zd a = 0, which means one may need a 
different frame for the energy-momentum and angular momentum components. For 
the purposes of obtaining directly a physical value for the observable quantities, one 
must separate the frame gauge dependent term into a physical energy plus a gauge 
dependent unphysical energy. This issue was first considered by Hechtj^^land he 
discovered a suitable remedy. The frame gauge dependent part can be separated 
using the identity 

£ z d a = di z d a + izdd a = Di z d a + i z Dd a - i z T a ^. (236) 

With the aid of this relation one can get frame gauge independent boundary terms 
for the quasi-local quantities. Thus, to drop the contribution from the frame gauge, 
one should make the replacement 

izT°p = D 0 Z a + {i z T a )p - {£ z d a )p DpZ a , (237) 

where we have introduced the convenient transposed connection: 

DZ a := DZ a + i z T a , (238) 

which naturally shows up whenever one expresses Lie derivative expressions in terms 
of a covariant derivative. In addition to our argument above, this replacement 
has been justified using (i) theoretical analysisj 49 ! 50 ! (ii) calculations for exact so¬ 
lutions,^ (iii) holonomic variables (for GR)P^land (iv) via a manifestly covariant 
Hamiltonian formulation.®® In our presentation here we used the coframe both for 
convenience and for its local translational gauge role; the reference just cited pro¬ 
vides a completely frame independent—manifestly covariant—alternative approach 
to the covariant Hamiltonian formalism. 

19.4. Quasi-local boundary terms 

With the aforementioned replacement we get our general set of symplectic quasi¬ 
local quantity boundary terms for the PG: 

B{Z) = i z 11“ } Ar a + Ad a A ** { !“ j 

+ {v»} A ^ + Ar “ eAiz {^} ■ (239) 

where the upper or lower line in each bracket is to be selected. As in thermody¬ 
namics, there are several kinds of “energy”, each corresponds to the work done in 
a different (ideal) physical process JM1I21 
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19.5. A preferred choice 

The cases of the PG that have been studied are mostly those where the first order 
potential is quadratic in the momentum fields, this leads to a linear relation be¬ 
tween the momenta and field strengths, and corresponds to 2nd order quasi-linear 
equations for the geometric potentials. The natural reference choice is Minkowsi 
spacetime, for which r M vanishes and D = D. For this class of theories, other things 
being equal, we would favor from among the set (123911 the Hamiltonian boundary 
term quasi-local choice (upper, lower, lower, upper), i.e., 

B{Z) = i z d a r a + AT a 0 A i zP J + DpZ a ApJ, (240) 

which leads to the following boundary term in the variation of the Hamiltonian: 

C(Z) = i z (6'd a At„- AT^pASpJ). (241) 

This corresponds to imposing boundary conditions on the coframe and the momen¬ 
tum conjugate to the connection. The associated energy flux expression is 

£zH(Z) = di z {£ z d a A t„ - AT a p A £ z pJ). (242) 

Regarding the total energy-momentum and angular momentum/center of mass 
momentum, the expression (124011 matches the expression (with Minkowski reference) 
for the PG Hamiltonian boundary term at spatial infinity that was first proposed 
by Hecht in 1993P^ 


19.6. Einstein’s GR 

Within the PG context, the special case of GR can be reached by imposing vanish¬ 
ing torsion with a Lagrange multiplier. In the general the first order formulation 
it is sufficient to take the potential to be independent of the coframe conjugate 
momentum r^. 

A first-order Lagrangian for vacuuir@ GR is 


£ G r = R af> A p a p + Dr A r M - V a0 A 


Pa/3 


2 K 


haft 


(243) 


which uses a Lagrange multiplier field V a ^ = to give the connection’s conju¬ 
gate momentum held a value which depends on the orthonormal frame: 


Pap y^Vap — 0 . 


(244) 


x For our purposes concerning the Hamiltonian boundary term, we consider here for simplicity just 
the vacuum case. The results will apply to all situations where the boundary of the region of 
interest is in the vacuum region, outside of the domain of the matter fields. That should include 
most of the cases of physical interest. 
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Variation of (I243D w.r.t. the coframe, connection and their respective momenta 
fields gives the (vacuum) equations: 

: 0 = Dt fj, + V ap A JL**,, (245) 

ST^ : 0 = Dp aP + d vl AT a] , (246) 

6t„ : 0 = D^, (247) 

Sp at 3 : 0 = R a/3 - V a/3 . (248) 

As expected (124711 gives vanishing torsion. From the differential of (124411 one gets 

Dp a/ 3 = A ri a f (249) 

Zk 

which vanishes, hence (12461) reduces to A r a i = 0, from which it follows that 
vanishes. Then (12451) with (12481) reduces to the vanishing of the Einstein 3-form: 

0 = -Ljrf* A Tja/Hfj, = -~G%p v . (250) 

By the way, had we included a suitable source, we would have obtained 

0 = ^R a p A p a \ + (251) 

where is the Hilbert energy-momentum 3-form. Using Dp a ^ ^ = 0 this relation 
can be rearranged as follows 

° = ^ [40 A + T a f, A d Va % + r a 7 A A 

- h A ~ r “ 7 A r7 ^ A A r7 ^ A ’*A] + v ( 252 ) 

The rearrangement identihes a certain superpotential 2-form, 

U M = ^r^A7 ?Q /3 /i , (253) 

and gravitational energy-momentum pseudotensor 3-form, 

2/st M = -r a 7 A F 7 ^ A rye' 3 ,, + A T 7 p A nJv (254) 

These manipulations and the resultant expressions are meaningful in both orthonor¬ 
mal and holonomic frames. In orthonormal frames they give the expressions for 
the so-called tetrad-teleparallel gauge current^^ while in holonomic frames we have 
obtained neat form versions of the Freud superpotential (1221) and the Einstein pseu¬ 
dotensor CoD. There is a remarkable contrast between the simple short calculation 
given here for these quantities and the long complicated ones discussed in section 0 
that were done in the past using tensor calculus. Via rearrangements of the held 
equations analogous to (12521) . generalized pseudotensors can be found for the PG.^ 
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19.7. Preferred boundary term for GR 

Over twenty years ago using the covariant Hamiltonian symplectic boundary term 
approach we proposed a quasi-local boundary term for GlP^ (an equivalent quasi¬ 
local expression obtained from a Noether argument using a global background ref¬ 
erence was proposed at about the same time by Katz, Bicak & Lynden-Be l 66 * 80 * ) : 

B(Z) = ^-(AT a 0 AizyJ+ D p Z a Ar 1 Jy, r?*" := *(0“ A & A • • •). (255) 

(This has the form of Hecht’s PG expression restricted to GR and natural extended 
to a boundary that need not be at infinity.) The boundary term in the variation of 
the Hamiltonian has the form 

5'H(Z) ~ di z (Ar a p A 5rj a P), (256) 

Since = *(d a AflP), this corresponds to fixing the orthonormal coframe (and 
thus the metric) on the boundary. The energy flux expression is 

£ z n(Z) = di z (AT a p A £ z yJ). (257) 

Like many other boundary term choices, at spatial infinity it gives the ADM,® 
MTWj^l Regge-Teitelboimf 1()7 Beig-0 Murchadhapl SzabadotP^J energy, momen¬ 
tum, angular-momentum, center-of-mass momentum. 

It has some special virtues: 

(i) at null infinity: the Bondi-Trautman energy & the Bondi energy fiuxpSl 

(ii) it is “covariant”, 

(iii) it is positive—at least for spherical solution^ and large spheres, 

(iv) for small spheres it is a positive multiple of the Bel-Robinson tensor,^ 

(v) first law of thermodynamics for black holesp^l 

(vi) for spherical solutions it has the hoop property,^ 

(vii) for a suitable choice of reference it vanishes for Minkowski space. 

20. A “best matched” reference 

In this section we turn to the second ambiguity that is inherent in quasi-local energy- 
momentum expressions: the choice of reference. Minkowski space is the natural 
choice, but one still needs to choose a specific Minkowski space. Recently we proposed 
(i) )D isometric matching on the boundary and (ii) energy optimization as criteria 
for selecting the “best matched” reference on the boundary of the quasi-local region. 


Note: for all other fields it is appropriate to choose vanishing reference values 
as the reference ground state—the vacuum. But for geometric gravity the standard 
ground state is the non-vanishing Minkowski metric, so a non-trivial reference is 
essential. One still needs to specify exactly which Minkowski space. 
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Reference values can be determined by choosing, in a neighborhood of the desired 
spacelike boundary 2-surface S, 4 smooth functions y 1 = y l (x 11 ), i = 0,1,2,3 with 
dy° A dy 1 A dy 2 A dy 3 ^ 0 and then defining a Minkowski reference by 

9 = -( dy°) 2 + (dy 1 ) 2 + (dy 2 ) 2 + (dy 3 ) 2 . (258) 

Geometrically this is equivalent to finding a diffeomorphism for a neighborhood 
of the 2-surface into Minkowski space. The associated reference connection is the 
pullback of the flat Minkowski connection: 

f a p = x a i (T i j yip + dy'p) = x a idy i p. (259) 

Here x a i is the inverse of y l p, where dy 1 = y l pd^. 

With these standard Minkowski coordinates y l , a Killing field of the reference 
has the form Z k = Zq + X^iy 1 , where A(f = Aq^, with Z§ and X kl being constants. 
The 2-surface integral of any one of our Hamiltonian boundary terms then has a 
value linear in these constant values: 

<f B(Z) = -Z k p k (S) + ^Ag%(S). (260) 

J s 

This implicity determines not only a quasi-local energy-momentum but also a quasi¬ 
local angular momentum/center-of-mass momentum. When specialized to GR the 
integrals p k (S), J k i(S) in the spatial asymptotic limit agree with accepted expres¬ 
sions for these quantities: MTW® §20.2 and Regge-Teitelboim 107 with the refine¬ 
ments of Beig-0 MurchadhgP and SzabadosP^ For the PG at spatial infinity, 
HechlPSI compared in detail his expression with the other proposed expressions, e.g., 
Ref. [ID]. At spatial infinity, with the asymptotics (11221) . all of our PG symplectic 
boundary terms (12391) will give the same values as Hecht’s expression 12401) . 

For energy-momentum one takes Z to be a translational Killing field of the 
Minkowski reference. Then the second term in our preferred quasi-local boundary 
expressions (124012551) vanishes @ With Z k = Z k = constant our preferred quasi-local 
expressions now take the form 

B pg (Z) = Z k x» k [T„ + (T a p - x a j dyip) A i^pj ], (261) 

B gk (Z) = Z k x^ k (T a p - x a j dy-’p) A . (262) 

20.1. The choice of reference 

To be completed, our Hamiltonian boundary term and the quasi-local energy- 
momentum/angular momentum proposal needs a prescription for choosing a ref¬ 
erence on the boundary. There are several options; one needs a choice suited to the 
application. 

For an extended region one may want a global background (see Refs. [10211031 
for this approach in GR). Consider for example solar system applications, more 

y For GR the second term in 12551 also vanishes for 4D isometric matching on S', a condition we 
shall use below. 
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specifically say we want to calculate using our quasi-local energy flux formula the 
tidal energy flux between Jupiter and its moon Io, that is believed to be responsible 
for Io’s volcanos. (This has already been done by several methods P^EHmil) On the 
other hand, if one wishes to study a given a metric expressed in a specific coordinate 
system the analytic approaclf 7 ^ may be a good choice. 

To explicitly determine the specific values of the quasi-local quantities one needs 
some good way to choose the reference. Minkowski spacetime is the natural choice, 
especially for asymptotically flat spacetimes. However, as noted above, almost any 
four functions will determine some Minkowski reference. With such freedom one 
can still get almost any value for the quasi-local quantities. This freedom is the 
quasi-local version of the second type of ambiguity mentioned in the introduction. 

Recently we proposed a progranPSI to fix the “best” choice for a quasi-local 
reference, i.e., one that is determined by the dynamical fields on the boundary. It 
has two parts: 4D isometric matching and optimization of a certain quantity. Here 
we present it along with a promising alternative optimization. 2 ^ For GR we have 
already found that our new procedure works well for an important special case: a 
certain class of axisymmetric spacetimes 127 which includes the Kerr metric. 

For the PG, not so much has been done yet. While PG energy-momentum and 
angular momentum calculations at both spatial and future null infinity were done 
long ago j 51 l 52 l and gave, in particular, the expected results for the Kerr metric, as far 
as we know there have not yet been any genuine quasi-local (i.e., for a finite region) 
PG calculations. This is not so surprising. In general the big obstacle is the 2D 
isometric embedding, which we are about to discuss. For spherical symmetry, all the 
calculations at least for GR can easily be done analytically. For the aforementioned 
class of axisymmetric metrics the 2D embedding problem has an algebraic solution. 
But the boundary expressions are already sufficiently complicated that the quasi¬ 
local energy integral for GR could only be evaluated numerically. Now that it is 
known how to do the case of axisymmetric GR the way is open for truly quasi¬ 
local PG energy and angular momentum calculations. For the PG, it seems that 
numerical calculations will be unavoidable. 

20.2. Isometric matching of the 2-surface 

We first recall an important procedure that has been used: isometric matching of 
the 2-surface S. This can be expressed in terms of quasi-spherical foliation adapted 
coordinates t, r, 9, ip as 

9ab = gAB = gijVAVB = -VaVb + SijVAVB , (263) 

where S is given by constant values of t, r, and A, B range over 9, ip. We use = to 
indicate a relation which holds only on the 2-surface S. Eq. (I263D is 3 conditions 
on the 4 functions y 1 . One can regard y° as the free choice. From a classic closed 
2-surface into R 3 embedding theorem—as long as S and y°( a^) are such that on S 

9ab := 9 ab + VaUb (264) 
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is convex—one has a unique embedding. Wang and Yau have discussed in detail 
this type of embedding of a 2-surface into Minkowski controlled by one function in 
their recent quasi-local work p 6 l 143 n 4 l 

Unfortunately, although there is a unique embedding, there is generally no ex¬ 
plicit analytic formula except in special simple cases, such as spherical symmetry or 
axisymmetry. The lack of an explicit formula for the solution of this 2D isometric 
embedding prevents exact quasi-local calculations in general. 

20.3. Complete 40 isometric matching 

Our “new” proposal is: complete 4-dimensional isometric matching on S: g=g : a 
part of which is still the just discussed 2D isometric embedding. 

In view of isometric matching, there should be a Lorentz transformation which 
on the boundary brings the dynamical coframe into line with the reference frame: 


if := L* a 0 Q =dy’. 

(265) 

The integrability condition for this equation is 

di? i |s = 0. 

(266) 


This 2-form equation restricted to the 2-surface gives 4 conditions on the 6 param¬ 
eter set of Lorentz transformations L l a . Thus 4D isometric matching has 6 — 4 = 2 
degrees of freedom. They can be identified as the choice of time embedding function 
y° in (12631) plus a boost parameter a in the plane normal to S. 

20.4. Optimal energy 

To fix the remaining 2 isometric matching parameters, one can regard the quasi¬ 
local value as a measure of the difference between the dynamical and the reference 
boundary values. This value will be a functional of y°, a. The critical points of this 
functional determine the distinguished choices for these 2 functions. 

Previously we proposed^ taking the optimal “best matched” embedding as the 
one which gives an extreme value to the associated invariant mass m 2 = —piPjg l C 
This should determine the reference up to a Poincare transformation. 

This is a reasonable condition, but, unfortunately, not so practical. The invariant 
mass is a sum of 4 terms, each quadratic in an integral over S. Note, however, that 
using the Poincare freedom one can get the same m value in the center-of-momentum 
frame from po- This leads us to our new proposal: take the preferred reference as 
one that gives a critical value to the quasi-local energy given by (12601) and (126211 or 
(12611) with Z k = Zq = <5q . We expect this much simpler optimization to give the 
same reference geometry as that obtained from using m 2 . 

Based on some physical and practical computational arguments it seems reason¬ 
able to expect a unique solution in general. In a numerical calculation in principle 

z For GR this was prop osed by Szabados at a workshop in Hsinchu, Taiwan in 2000. He has since 
investigated it in detail 
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one could just calculate the energy values given by (12601) and (12611) or (I262D with 
Z k = i5q for a great many choices of y°,a subject to the 4D isometric matching 
constraint (12651) and the integrability condition (12661) . and then note the energy 
critical points. 

For GR this “best matching” procedure already gave reasonable quasi-local en¬ 
ergy results for spherically symmetric systems j 2 1 l 79 l 149 l U Hl p or the Schwarzschild 
metric the “best matched'’ quasi-local energy has the value first found by Brown 
and YorkjSl 


E(r) 


2 m 

1 + %JT^2rnJr' 


(267) 


using a closely related boundary term and a 2-surface embedding into R 3 . Now 
we also have sensible results for certain axisymmetric systems including the Kerr 
metricFor the surface of constant Boyer-Linquist r, the angular momentum is 
simply a constant independent of r, equal to its usual asymptotic value, J. The 
quasi-local energy integral, however, is not so simple and can only be evaluated 
numerically. 


21. Concluding discussion 

In addition to her two key theorems regarding symmetry, Emmy Noether in her 1918 
paper also proved that for diffeomorphically invariant gravity there is no proper to¬ 
tal energy-momentum density. In other words there is no covariant total energy 
momentum density tensor for gravitating systems. In Ch. 20 of their textP^ Mis- 
ner, Thorne and Wheeler discuss this feature as a consequence of the equivalence 
principle and argue that only the total energy-momentum of gravitating systems 
is meaningful. But clearly the gravitational interaction is local and does allow for 
the local exchange of energy-momentum. To account for this various non-covariant 
expressions called pseudotensors (each generated by a certain superpotential) have 
been proposed. There thus arose two ambiguities: which expression? in which 
reference frame? The modern idea is quasi-local: energy-momentum is associated 
with a closed 2-surface. 

One approach, which is the one we have used, is via the Hamiltonian. With the 
aid of differential forms and a first order variational formulation, we have developed 
a covariant Hamiltonian formalism. The Hamiltonian boundary term plays key 
roles: it determines the boundary conditions and the quasi-local values. We have 
shown that this approach includes all the traditional pseudotensors while clarifying 
the ambiguities: each superpotential is a possible Hamiltonian boundary term which 
is associated with a specific boundary condition, and the reference frame becomes 
a choice of the reference values (ground state) on the boundary. 

One can regard gravity as a local gauge theory for the global symmetry group for 
Minkowski spacetime: the Poincare group. The appropriate geometry is Riemann- 
Cartan: the curvature is the field strength for Lorentz transformations and the 
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torsion is the field strength for local translations (infinitesimal diffeomporhisms). 
For comparison, we included in our discussion a general internal gauge gauge field. 
In this way we can identify the analog of the internal gauge vector potential. For 
the local Lorentz symmetry it is the (metric compatible) connection one-form', for 
the local translation symmetry it is the orthonormal coframe. We developed in 
considerable detail the associated Lagrangian and Hamiltonian Noether currents 
and differential identities, so one can compare the similarities and differences of the 
spacetime gauge symmetry expressions with those of a generic internal symmetry. 
Briefly, the Lorentz rotation sector is quite similar to that of an internal gauge sym¬ 
metry, but the translation symmetry has both striking similarities and differences. 
This would be much less clear if we had opted for a formulation which does not 
include the coframe as a dynamical variable. In the approach of Blagojevic & HehlP 
the use of the orthonormal frame is motivated by the need to describe spin; here 
our basic motivation is in terms of the coframe’s fundamental gauge role regarding 
local spacetime translations. 

The geometric/gauge symmetry approach is helpful in identifying a good ex¬ 
pression for our Hamiltonian boundary term expression for quasi-local quantities. 
Our preferred expression for GR turns out to correspond to fixing the metric on the 
boundary—which is obviously a reasonable boundary condition choice. 

A notable feature of the Hamiltonian boundary term for dynamic geometry is 
that it necessarily depends on the choice of some non-dynamical reference values 
(this is a manifestation of Noether’s result regarding non-existence of a proper 
energy-momentum density). One can get almost any quasi-local value if one allows 
free rein in the choice of reference. One needs to fix a non-dynamical reference frame, 
but only on the boundary of the region. This can be compared to choosing some flat 
plane to map a part of the curved surface of the Earth. One could slice the plane 
through the surface of the Earth; for a spherical Earth the planar geometry would 
exactly match on a circle. Similarly for spacetime. On the 2-dimensional boundary 
of a spatial region one can exactly match the curved 4D spacetime metric with a flat 
Minkowski spacetime metric. Detailed analysis shows that there is still two degrees 
of freedom. We proposed that a good way to fix these was by considering the critical 
points of our quasi-local expression. There might be other sensible options, but the 
main point is that a reasonable choice is available. 

Our principal results, the Hamiltonian boundary terms that are our preferred 
quasi-local energy-momentum expressions (1240I255[) for the PG and GR, were ob¬ 
tained by considering the Hamiltonian, geometry, Noether symmetry, and space- 
time gauge theory. The harmonious combination of all of these perspectives makes 
a strong case for the results. Nevertheless, it should be noted that one can also be 
led to this result from other perspectives. Regarding GR, essentially the same ex¬ 
pression has been obtained (i) via a Noether approach with a global reference) 66 * 80 ! 
and (ii) via a symplectic covariant Hamiltonian approach using the metric in a 
holonomic frameFor the PG (including GR as a special case), the same pre- 
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ferred expression was found via an entirely frame independent manifestly covariant 
Hamiltonian formalism.^ Although in principle there are an unlimited number of 
possible Hamiltonian boundary term quasi-local energy-momentum expressions— 
which are in a formal sense all of equal status- in practice one can—with very good 
reasons—discover that one expression stands out. 
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